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MODEL FOR MICROFILTRATION OF POLY-DISPERSE 

SUSPENSIONS 

[00011 This application claims the benefit of U.S. Provisional Patent 

Application Serial No. 60/403,575, filed August 14, 2002, and U.S. Provisional 
Patent AppUcation Serial No. 60/471,603, filed May 19, 2003. 



FIELD OF THE INVENTION 

[0002] The present invention relates to a model for microfiltration of poly- 

10 disperse suspensions and solutions. 

BACKGROUND OF THE INVENTION 

[0003] There are presently maay solutions containing macromolecular and 

suspended particles tihat are being filtered with synthetic membranes. These 

15 include waste water, siuface water, automotive paint streams, and a host of 
streams firom the biomedical and bioprocessing industries containing proteins, 
cells, DNA, fat globviles, colloids, milk, suspended particles etc,. Microfiltration 
of cell culture, fermentation broths, blood, and milk are representative examples 
of streams firom the latter industries. (Belfort et al., *The Behavior of Suspensions 

20 and Macromolecular Solutions in Crossflow Microfiltration," J. Membr. Scu 96:1- 
58 (1994); Nagata et al., "Qrossflow Membrane Microfiltration of a Bacterial 
Fermentation Broth," BiotecftnoL Bioeng. 34:447-466 (1988)). Because of the 
complex nature of these fluids, and due to the difficulty in specifying the behavior 
of their suspended and dissolved components, modeling miarofiltration has been 

25 difficult. Belfort and coworkers have summarized the behavior of suspended 
particles during microfiltration. (Green et al., 'Touling of Ultrafiltration 
Membranes: Lateral Migration and the Particle Trajectory Model," Desalination, 
35:129-147 (1980); Altena et al., '"Lateral Migration of Spherical Particles in 
Porous Channels: Application to Membrane Filtration," Chem. Eng. Sci. 39:343- 

30 355 (1984); Alt^ia et al., "Lateral Migration of Spherical Particles in Porous Tube 
Flows: Channels: Application to Membrane Filtration," Physico-Chem. Hydrodyn 
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6:393-413 (1985)). Oflier researchers have studied the behavior of 
macromolecules diuing microfiltration (Zydaey et al., "A Concentration 
Polarization Model for the Filtrate Permeation Flux in Crossflow Nficrofiltration 
of Particulate Suspensions," Chem, Eng. Commun. 47:1-21 (1986); Rautenbach et 
5 al., ^Ultrafiltration of Macronaolecxilar Solutions and Crossflow Microfiltration of 
Colloidal Suspensions. A Contribution to Permeate Permeation Flux 
Calculations," J. Membr. Set 36:231-242 (1988); Samuelsson et al., *Tredicting 
Limiting PCTneation Flux of Skim Milk in Crossflow Microfiltration," J. Membr. 
Sci. 129:277-281 (1997); Romero et al., "Global Model of Crossflow 

10 Microfiltration Based on Hydrodynamic Particle DiflEusion," /. Membr, Sci. 

39:157-185 (1988)). Only a few studies have been reported in the literature on 
modeling the behavior of poly-disp^e feeds containing both macromolecules and 
suspended particles for microfiltration (Dharmappa et al., "A Comprehensive 
Model for Crossflow Filtration Incorporating Polydispersity of the Influent," J. 

15 MembK ScL 65:173-185 (1992); Ould-Dris et al., ^TEffect of Cake Thickness and 
Particle Polydispersity on Prediction of Permeate Permeation Flux in 
Microfiltration of Particulate Suspensions by a Hydrodynamic DifiRision Model," 
J. Membr. ScL 164:211-227 (2000)). 

[0004] Crossflow (also known as tangential) microfiltration is a very 

20 complex phenomenon. Some of tixe variables which influence permeation flux 
and retention are membrane type and chenustry, module geometry, particle size 
distribution, nature of particles, interaction between individual particles and 
between particles and the membrane, fluid dynamics, operating mode, 
temperature, and pH and ionic strength of the media. This is a formidable set of 
25 variables and, to date, no unified theory exists to give a rigorous expression of 
permeation flux and retention for crossflow microfiltration. Existing theories 
render the problem tractable by concentrating on, at best, only a few aspects of the 
problem. Therefore, for any given case, different flieories may yield widely 
divergent results. In applying a theoretical model, extreme care must be exercised 
30 to check the specifics of the case and critically evaluate the dominant parameters 
and compare these with the model assimiptions. It is quite possible, particularly in 
dealing with complex suspensions, that no one phenomenon or class of 
phenomena is dontiinant In such a case, an existing model may not give the true 
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picture and a new model may need to be evolved to capture ihe dominant 
phenomena. Usually, with the onset of microfiltration at constant transmembrane 
pressvire, there is a rapid decline of permeation flux due to concentration 
polarization and pore constriction, followed by a quasi steady state where there is 
5 a gradual decline in permeation flux due to particle deposition and increase in 
paiticle concentration (and hence viscosity) of the bulk solution (Belfort et al., 
*The Behavior of Suspensions and Macromolecular Solutions in Crossflow 
Mictofiltration,'V. Membr. Sci. 96:1-58 (1994); Nagata et al,, "Crossflow 
Membrane Microfiltration of a Bacterial Fermentation Broth," BiotechnoL 

10 Bioeng, 34:447-466 (1988)). Several models exist which attempt to evaluate the 
quasi steady state permeation flvix. These models are based on equilibrium 
between transport to and back-transport of particles firom the membrane wall. It is 
important to note that these models are valid in flie laminar flow regime for fully 
retentive mraibranes, deal with the pressure-independent permeation flvix regime, 

1 5 and ignore particle interaction with the membrane. Prediction of permeation flux 
in the pressure-dependent and transient regimes has been addressed by Romero 
and Davis (Romero et al., "Global Model of Crossflow Microfiltration Based on 
Hydrodynamic Particle Diffusion," J. Membr, ScL 39:157-185 (1988); Romero et 
al., ^Transient Model of Crossflow Microfiltration," Chem. Eng. ScL 45:13-25 

20 (1990)). A simplified version of the gel-polarization model for describing 
concentration polarization in ultrafiltration for fully retentive membranes, 
originally presented by Blatt et al,, "Solute Polarization and Cake Formation in 
Membrane Ultrafiltration: Causes, Consequences and Control Techniques," in 
Membrane Science and Technology^ J. E. Flinn, ed.. New YorkrPlenum, pp. 47-97 

25 (1970), is given by the permeation flux at an axial position x along flie flow path, 

J(x) = /fc(x)ln(<2i^/(^b) (1) 

where the membrane and cake resistances do not appear in this formxilation and 
30 k(x), the mass transfer coefficient is given by the ratio of the molecular diffusion 
coefficient, D, to the x - dependent mass boimdary layer thickness, <5(x). 
[0005] Using the Leveque solution for laminar flow in a closed tube and 

the Stokes-Einstein relationship for Brownian diffusion of mono-dispersed 
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Spheres of radius a, Eq. (1) changes to Eq. (2), shown in Table 1. The Brownian 
(molecular) diflEusion model is valid for particle diameters below --0.1 jxm and at 
low axial shear rates. For microfiltration situations where a> I |jm, Eq. (2) 
under-predicts the permeation flux by one to two orders of magnitude. Green and 
5 Belfort (Green et al., *Touling of Ultrafiltration Membranes: Lateral Migration 
and the Particle Trajectory Model/' Desalination^ 35:129-147 (1980)), referring to 
this as the flux paradox, searched for mechanisms that induce back-migration 
from the membrane surface in addition to molecular diffusion. Belfort and 
coworkers invoked inerdal lift as an additional back (or lateral) migration 

10 mechanism and used it to explain the disorepancy in permeation flux mentioned 
above (Green et al., 'Touling of Ultrafiltration Membranes: Lateral Migration and 
the Particle Trajectory Model,*' ZJesa/z/wtfo/i, 35:129-147 (1980); Altena et al., 
'*Lateral ^/figration of Spherical Particles in Porous Channels: Application to 
Membrane Filtration,*' Chem, Eng, Set 39:343-355 (1984); Altena et al., "Lateral 

15 Migration of Spherical Particles in Porous Tube Flows: Chaimels: Application to 
Membrane Filtration," PAy5zcc?-cAem. Hydrodyn. 6:393-413 (1985)). Their 
expression for permeation flux is given by Eq. (3) in Table 1. This model is 
applicable for dilute solutions containing large particles (> 20 |im) and high axial 
shear stress or fast laminar flow situations. Zydney and Colton (Zydney et al., "A 

20 Concentration Polarization Model for the Filtrate Permeation Flux in Crossflow 
Microfiltration of Particulate Suspensions," Chem. Eng. Commun. 47:1-21 (1986)) 
resolved the permeation flux paradox by proposing that the gel-concentration 
polarization model should be used, provided the molecular diffusion term in the 
Leveque solution is replaced by the shear-induced hydrodynamic diffusivity (and 

25 ttie permeation flux is averaged over the axial flow path), first measured by 
Eckstein et al, "Self-Diffusion of Particles in Shear Flow of a Suspension,'' J. 
FluidMech. 79:191-208 (1977)). This model, given by Eq. (4) in Table 1, has 
been experimentally validated for latex and blood suspensions for a broad range of 
conditions. It is suitable for particles in the diameter size range 0,1 to 20 pm. Li 

30 et al., "An Assessment of Depolarization Models of Crossflow Microfiltration by 
Direct Observation tiurough flie Membrane," J, Membr, Set 172:135-147 (2000)), 
based on direct observation of particle dynamics during membrane filtration. 
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suggested that the coefficient in Eq. (4) should be replaced with 0.0595. Using a 
similarity solution, Davis and Sherwood (Davis et al., "A Similarity Solution for 
Steady-State Crossflow Microfiltration," Chem. Eng. ScL 45:3203-3209 (1990); 
Ho et al., Membrane Handbook, New York: Van Nostrand Reinhold; pp 480-505 
5 (1 992)) have accounted for a concentration-dependent flxiid viscosity and 
difJEiisivity and obtained Eq, (5), shown in Table 1 . The results from this 
e3q)ression agree with those from Eq. (4) to within 15% for 0.01 and 0.6, 
which is close to the maximum for monodisperse rigid spherical particles. 
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[0006] The models discussed above deal with crossflow microfiltration of 

idealized mono-disperse streams, whereas most real microfiltration applications 
deal with poly-disperse feed streams with particle sizes ranging iBrom the 
macromolecular to the colloidal. This has generated the need to develop a 
generally appUcable model useful for predicting length-averaged permeation flux 
and yield of target molecules for poly-disperse feeds imder a variety of filtration 
conditions. 

[0007] The present invention is directed to overcoming these and other 

deficiencies in ihe art. 



SUMMARY OF THR TNVENTION 

[0008] The present invention relates to a method for predicting pressure 

independent permeation flux and target molecule yield in a permeate resulting 
firom crossflow membrane filtration of particles in a poly-disperse suspension. 
This method involves determining the particle size distribution of the poly- 
disperse suspension, the equivalmt spherical raxiii of the particles, the viscosity of 
the suspension, and the maximum back-transport velocity (wj) for all particles. It 
also involves estimating the maximimi aggregate packing volume firaction {^m) for 
all particles at a wall of the filtration membrane firom geometric considerations, 
and selecting the particle that gives a minimum permeation flux at a given 
filtration membrane shear rate, where the selected particle has a radius (aO, and 
determining a predicted permeation flux. The mefliod also involves determining 
packing density at ^ a membrane wall for each particle size (aj for J ^ z ) at the 
predicted permeation flux. Also determined are interstitial packing dmsity 
i^iinterstice) of particlcs in the suspension which are smallest, and minimum pore 
diameter (2rminimum) based on the packing density of each particle. The yield of a 
target species in the filtration permeate is then estimated by calculating observed 
sieving coefficient (So) for the target species. As a result, permeation flux and 
target molecule yield of the poly-disperse suspension during crossflow filtration 
are predicted. 
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[0009] The present invention also relates to a method for detennining the 

packing density of particles of a poly-disperse suspension at a naembrane wall. 
Tbis method involves providing a predicted permeation flux (/), dete rmin i n g the 
packing density for all particle sizes at the predicted permeation flux, and 
detennining interstitial packing density {^interstice) of particles in the suspension 
which are smallest, thereby detennining the packing density at the membrane wall 
of particles of the poly-disperse suspaision. 

[0010] Another aspect of the present invention is a method for predicting 

pressure independent permeation flux for crossflow membrane filtration of a poly- 
disperse susp^ion* This method involves determining the viscosity of the 
suspension, detennining the maximum back-transport velocity (ui ) for all 
particles, and estimating the maximum aggregate packing volume fraction 
all particles at a wall of the filtration membrane wall from geometric 
considerations. The particle that gives a minimirai permeation flux at a given 
filtration membrane shear rate is selected, where the selected particle has a radius 
(aO- A predicted permeation flux (J) is determined, and the packing density (^J) 
at the membrane wall for each particle size (o/ for j ^ i) is determined at the 
predicted permeation flux. 

[0011] Another aspect of the present invention is a method for calculating 

yield of a target molecule in a permeate for a poly-disperse suspension during 
crossflow membrane filtration. This method involves determining the minimum 
pore diameter (2rnuniniuni) based on the packing density of each particle, and 
estimating the yield of a target species in the filtration permeate by calculating the 
observed sieving coefficient (So) for the target species. 
[0012] The present invention also relates to a method for designing a 

crossflow membrane filtration system for a poly-disperse suspension. This 
method involves selecting a poly-disperse suspension and predicting pressure 
independent permeation flux and target molecule yield in a permeate resulting 
from crossflow membrane filtration of particles in a poly-disperse suspension as 
described above. Conditions for filtration based on tihie prediction of permeation 
flux and target molecule yield are then optimized to design a filtration system for 
the selected poly-disperse suspension. 
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[0013] Yet another aspect of the preset invCTition is a method of selecting 

operating conditions of a crossflow filtration system for poly-disperse 
suspensions. This method involves predicting the pressure independent 
penneation flux and target molecule yield in a permeate resulting from crossflow 
membrane filtration of particles in a poly-disperse suspension as described above. 
Operating conditions of the system are selected using the limiting pressure 
independent penneation flux determined for a given shear rate to obtain an 
optimal balance between permeation flux and yield of a target species. 
[0014] The present invention also relates to a method of modeling a 

process for filtration of a poly-disperse suspension. This method involves 
applying the method for predicting pressure independent permeation flux and 
target molecule yield in a permeate resulting firom crossflow membrane filtration 
of paiticles in a poly-dispCTse suspension, as described above, and using a 
computer-generated program to model a process for filtration of a poly-disperse 
suspension, 

[0015] Most microfiltration applications deal with poly-disperse feed 

streams with particle sizes ranging from the macromolecular to the colloidal. This 
has generated the need to develop a model which can predict the pressure 
independent permeation flux, a priori, and can provide an estimate of the yield of 
target molecules (i.e., the product concentration in the permeate) in the permeate 
for poly-disperse suspensions. Thus, the present invention relates to a method for 
calculating yield of a target molecule in a permeate for a poly-disperse sxispension 
during crossflow membrane filtration. The iterative methodology is capable of 
predicting crossflow microfiltration performance of poly-disperse suspensions in a 
variety of settings. For example, it is suitable for filtration of a poly-disperse 
suspension of waste water, surface water, paints, metallurgical wastes, 
environmental pollutants, industrial waste streams and streams from biochemical 
and bio-processing industries containing proteins, cells, DNA, colloids, milk, 
suspended paiticles etc. 

[0016] Several researchers have reported that there exists a critical ratio 

between permeation flux and wall shear rate beyond which protein transmission 
drops drastically and permeation flux does not increase with increasing 
transmembrane pressure (Berre et al., "Skim Milk Crossflow Microfiltration 
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Performance Versus Penneation Fliix to Wall Shear Stress Ratio," J. Membr. ScL 
117:261-270 (1996); Gesan-Guiziou et al., "Critical Stability Conditions in 
Crossflow Microfiltration of Skimmed Milk: Transition to Irreversible 
Deposition," J. Membr. ScL 158:21 1-222 (1999); Gesan-Guiziou et al., 
"Performance of Whey Crossflow Microfiltration During Transient and Stationary 
Operating Conditions,^'/. Membr. Set 104:271-281 (1995); Gesan-Guiziou et al., 
"Critical Stability Conditions in Skimmed Milk Crossflow Microfiltration: Impact 
on Operating Modes," Lait 80: 129-140 (2000)). The method presented here 
provides a quantitative firamewoik to explain this phenomenon and also provides 
tiie rationale for operating at uniform axial transmembrane pressure along tiie flow 
path. 

[00171 For simplicity, the adhesion of particles to the surface of the 

membrane, charge interactions between particles and the membrane, and between 
particles, are ignored. Several papers have demonstrated that electrostatic effects 
actually can increase permeation flux during filtration (Zydney et al., 
"Intermolecular Electrostatic Interactions and Their Effect on Permeation Flux 
and Protein Deposition During Protein Filtration," Biotechnol Prog. 10:207-213 
(1994); Aimar et al., "Fouling and Concentration Polarization in Ultrafiltration 
and Microfiltration," Membrane Processes in Separation and Purification, Crespo 
et al., eds., Dordrecht, Boston, London:Kluwer Academic Publishers (1994)). 
Under these conditions, this renders flie current approach conservative and 
compensates to some extent for neglecting particle aggregation and particle 
adhesion to the membrane. All particles are assumed to be spherical, with 
volumes equal to the actual particle volumes. 

[0018] The following scenario is proposed. There are three regimes of 

operation, as shown schematically in Figure 1 A. In the first regune, a rapid 
decline of permeation flux occurs with respect to the water flux, mainly due to 
deposits within the membrane pores leading to flow constriction. The deposits on 
the membrane have not yet consolidated. Instead, the deposits form a sparse layer 
which caimot be compressed. The permeation flux is proportional to 
transmembrane pressure but the slope of tiie Ime corresponds to the inverse of the 
resistance (or to the permeability coefficient) of the constricted and not the clean 
membrane. This results in a penneation flux that is lower than the clean water 
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permeation flux. This approach tiierefore deviates from that suggested by Howell 
et al. (Howell, "Sub-Critical Permeation Flux Operation of Microfiltration," J. 
Membr. ScL 107:165-171 (1994); Field et al., "Critical Permeation Flux Concept 
for Microfiltration Fouling," Membn ScL, 100:259-272 (1994)), who have 
proposed the existence of a critical permeation flux in microfiltration below which 
there is no fouling at all. With an increase in transmembrane pressure, the 
permeation flux along with the concentration of deposits increases until the gel 
(Landman et al., '^Pressure Filtration of Flocculated Suspensions," ^CAS J. 
41:1687-1700 (1995)) concentration is exceeded. The gel concentration is that 
threshold concentration beyond which the deposits consolidate into an 
intercounected network enable of becoming denser with increasing pressure. Jn 
this regime, an increase in transmembrane pressure results not only in an increase 
in the permeation flux but also an increase in the cake resistance due to cake 
drasificatLon. Tims, the permeation flux-transmembrane pressure curve falls 
below and away from the initial straight line. Li the third regime, the permeation 
flux does not increase with increasing transmembrane pressure. Layers of 
particles develop on the membrane surface concomitant with a deposition of 
particles within the pores. The particle layer adjacent to the membrane will 
correspond to the maximxnn packing density. This could range from 0.64 for rigid 
monodisperse spheres to 0.95 for highly deformable particles like red blood cells 
(Zeman et al., "Microfiltration and Ultrafiltration Principles and Applications, 
**New York:Marcel Dekker, Inc. (1996) or poly-disperse rigid sph^es of widely 
different sizes. The deposition within the pore depends on the membrane 
chemistry. For hydrophilic membranes particle adhesion will be lower than for 
hydrophobic membranes for most proteins (Hanemaijier, ^'Microfiltration in Whey 
Processing," Z)esa//watfo« 53:143-155 (1985); Koehler et al., 'Tntermolecular 
Forces Between Proteins and Polymer Films with Relevance to Filtration," 
Zangmmr 13(15):4162-4171 (1997)). This is because a hydrophilic membrane 
has a greater affinity for water in comparison with proteins, leading to the 
preferential formation of water layers, instead of protein layers, adjacent to the 
membrane. If the transmembrane pressure is increased, the particle layers above 
the compact layer adjacent to the membrane become denser. Therefore, the cake 
resistance increases and the permeation flux is maintained. This leads to a quasi 
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steady state permeation flux as long as the bulk suspension remains the same and 
the operating conditions are maintained. There is a gradual decline in permeation 
flux due to the increase in deposits within the pores leading to higher pore 
constriction and a time dependent compression of the cake layer. Essentially, the 
cake forms due to the balance of convective transport of solutes towards the 
membrane with back transport of solutes jBrom the membrane to the bulk and the 
solute transport through the membrane. The cake and the constricted membrane 
provide the necessary resistances to maintain the permeation flux required to 
achieve this balance. 

[0019] A theory of cake filtration relating to the composition and sieving 

characteristics of the deposited cake has been developed based on an experimental 
determination of a rheological parameter related to the specific cake resistance. 
(Landman et al., *Tressure Filtration of Flocculated SuspGnsions,^* AIChE J. 
41:1687-1700 (1995)). Other researchers have experimeatally evaluated the 
sieving characteristics of filter cakes in microfiltration formed by protein deposits 
(Palecek et al., **Hydraulic Permeability of Protein Deposits Formed During 
Microfiltration: Effect of Solution Ph and Ionic Strength," X Membr. Sci. 95:71- 
81 (1994); Meireles et al., "Effects of Protein Fouling on the Apparent Pore Size 
Distribution of Sieving Membranes/ V. Membr. ScL 56:13-28 (1991) (NOTE: 
This paper has a definitional error — for correction see Zydney et al., *TJse of the 
Log-Normal Probability Density Function to Analyze Membrane Pore Size 
Distributions: Functional Forms and Discrepancies," Membr, ScL 91:293-298 
(1994)); Mochizuki et al,, "Sieving Characteristics of Albumin Deposits Formed 
During Microfiltration," J. Colloid Interface Sci. 158:136-145(1993)). In this 
invention, the predominant nature of the cake is considered to be determined by 
the propensity of particles of different sizes to difRise back to the bulk for a given 
shear rate. The cake is thought to be primarily composed of particles with the 
least back-transport velocity. This does not preclude the possibility of smaller 
particles firom lodging themselves in the interstices of a compact layer of larger 
particles, as shown in Figm'e IB. The packing density and nature of the cake will 
correspond to that necessary, at equilibrium, to support a combination of back- 
transport rate and solute transport through the cake to match the convective 
transport of solutes to the membrane wall (i.e. the balance condition). Although 
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the packing density of the smaller paiticles could be less than that of the larger 
particles, the interstices are smaller and could very well determine the extent of 
passage of particles through the cake. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] Figures 1 A-C are schematic representations of various parameters 

that influence the microfiltration of poly-disperse suspensions. Figure 1 A is a 
schematic of three different operating regimes during microfiltration of poly- 
disperse suspensions. Regime I is pore constriction. Regime n is cake 
consolidation. Regime III is pressure-independent flux. Figure IB is a schematic 
of a sparse cake comprising a bi-disperse mixture of large and small particles 
while operating in Regime I. Figure IC shows particle packing constraints for the 
test case of a face-centered cubic arrangement (explained in greater detail in 
Example 3): ^ + <fc + ^ <|)Max ; <^ + ^ constant (0.68 in this case) and <fe < 
0.74[1 — (^ + jfc)]. Transmission of particle of radius a\^t.O for ^terstice ^ 0.74. 
[0021] Figures 2A-B are graphs of predicted values for a hypothetical 

suspension comprising particles of three sizes. Figure 2A is a graph of predicted 
values of pressure-independ^t permeation flux based on back-transport of 
particles of 10 imi (solid line), 180 nm (short dashed line), and 300 nm (long 
dashed line) versus mean axial shear rate. Figure 2B is a graph of predicted 
values of poly-disperse pressure-independent permeation flux versus mean axial 
shear rate, for a 6-fiber hollow fiber module of length 300 mm, internal diameter 
1 .27 mm, and pore diameter of 1 00 nm filtering whole transgenic goat milk at 298 
K. 

[0022] Figure 3 is a graph of predicted values of pressinre-independent 

permeation flux versus mean axial shear rate for 6-fiber hollow fiber modules of 
various lengths. The internal diameter and pore diameter of the fibers was 1.27 
mm and 100 nm, respectively. 

[0023] Figure 4 is a graph of predicted values of pressure-independent 

permeation flux versus mean axial shear rate for a 6-fiber hollow fiber module of 
length 300 mm, internal diameter 1 .27 mm, and pore diameter of 100 nm at 
various bulk feed volume firactions. 
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[0024] Figure 5 is a graph of predicted values of pressure-independeat 

permeate flux and yield of 10 mn particles versus mean axial shear rate after 
diafiltration of 4 diavolumes for a 6-fiber hollow fiber module of length 300 nmi, 
internal diameter 1.27 mm, and pore diameter of 100 nm. 
[0025] Figures 6A-C are a schematic of variation in pressure inside the 

bore of a hollow fiber (PI to P2) with varied length. Figure 6 A shows the 
permeate pressxjre in the extra-capillary space is constant at P3, wh^e beyond an 
axial distance L% the permeate pressure exceeds the retentate pressure causmg 
penneate to flow in the reverse direction firom the permeate side to the retentate 
side, called Starling flow. Figure 6B shows the constant at P3 for a short axial 
path length such that P3<P2. In this case reverse permeation does not occur 
because P2 > P3, as shown in Figure 6C, is allowed to decline firom PSA to P3B in 
the extra-capillary space as the pemxeate stream is pumped axiaUy in co-flow 
mode with respect to the retentate stream such that the transmembrane pressure is 
kept constant wifli axial distance. 

[0026] Figure 7 is a diagram of the linear and new coiled hollow fiber 

design according to U,S^ Patent No. RE 37,759 to Belfort, which is hereby 
incorporated by reference in its entirety. 

[0027] Figure 8 is a flow diagram of the dual hollow fiber microfiltration 

test system used for microfiltration of transgenic goat milk. 
[0028] Figures 9A-B are graphs of predicted pressure-independent 

permeation flux values. Figure 9A is a graph of predicted values of pressure- 
independent permeation flux based on back-transport of particles of IgG (solid 
line), casein micelles (short dashed line), and fet globules Gong dashed line) 
versus mean axial shear rate. Figure 9B shows pressure-independent poly- 
disperse permeation flux versus mean axial shear rate for a 6-fiber hollow fiber 
module of length 300 mm, internal diameter L27 rmn, and pore diameter of 100 
nm filtering whole transgenic goat milk at 298 K. 

[0029] Figure 10 is a graph of the friction factor versus Reynolds number 

for linear (fiill symbols) and helical (open symbols) with DI water. The solid line 
passing through the data for the helical modules is the Mishra and Gupta equation: 
fheiicai = (64/Re)[l+ 0.0033(logioDe)4] ; and the dashed Une passing tiurough the 
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data for the linear modxile is the Hagen- Poiseuille equation : linear = 64/Re. AL- 
1, ♦ L-2, ■ L-3, A H-1, 0 H-2, □ H-3. 

[0030] Figure 1 1 is a graph of permeation flux as a function of 

transmembrane pressure ("TMP") for the linear module (30 cm long) at different 
whole transgenic goat milk bulk protein concentrations (■ 13 g/1; ♦ 23 g/1; • 42 
gA; A 54 g/1; X 68 g/l). The axial Reynolds number was kept at 1500. The near 
vertical line (y = 37x) is the praneation flux with DI water, 
[0031] Figure 12 is a graph showing permeation flux as a function of 

transmembrane pressure ('TMP") for the helical module (24.9 cm long) at 
different whole transgenic goat milk bulk protein concentrations (0 13 g/1; 0 23 
g/l; o 42 g/1; A 54 g/1; x 68 g/1). The axial Reynolds number was kept at 1075. 
The near vertical line (y = 33x) is the permeation flux with DI water. 
[0032] Figure 1 3 is a graph of permeation flux versus the logarithm of the 

protein feed concentration for various whole transgenic goat nadlk microfiltration 
runs for the linear (♦) and the helical (0) modules. The plot is made according to 
the gel polarization model J = k Ln(Cw/Cb). The resulting linear equations are y 
= -26.5x + 161,8, r2 = 0.99 for the helical module indicating kheUcai = 26.5 knh, y = 
-18.2x + 1 10.9, r2 = 0.95 for the linear module indicating kunear = 1 8.2 Imh. Based 
on the x-intercept, Cw (helical) = 448 g/1 and Cw (linear) = 443 g/1. The 
transmembrane pressures corresponding to the pressure-independent flux were 
about 60 kPa for the linear and 70 kPa for the helical modules. 
[0033] Figure 14 is a graph of permeation flux versus number of 

diavolumes passed through the membrane for different starting concentrations of 
whole transgenic goat milk at 298 K and Reynolds number of 1400. The various 
cases are ■ Linear IX, • Linear 1.5X, ♦ Linear 2X, ALinear 3X, □ Helical IX, o 
HeUcal 1.5X, 0 Helical 2X, x HeUcal 2.5X, A HeUcal 3X. 

[0034] Figure 15 is a graph of IgG yield versus whole transgenic goat milk 

starting concentration (X) for helical (open symbols) and linear (full symbols) 

modules at 298 K and Reynolds number of 1400 after 5 diavolumes. 

[0035] Figure 16 is a graph of observed protein sieving coefficient versus 

bulk protein concentration of whole transgenic goat milk for helical (open 

symbols) and linear (ML symbols) modules at 298 K and Reynolds number of 

1400. 



wo 2004/016334 



-16- 



PCT/US2003/025230 



[0036] Figure 17 is a graph of permeation flux versus number of 

diavolumes passed trough the membrane for 2x concentration of whole 
transgenic goat mUk at 298 K and various Reynolds numbers. The various cases 
are ▲ Linear Re = 860, • Linear Re = 1285, ■ Linear Re = 1500, ♦ Linear Re = 
1715, X HeUcal Re = 643, A Helical Re =860, o HeUcal Re =1070, □ HeUcal Re 
=1400, 0 HeUcal Re =1610. 

[0037] Figure 18 is a graph of permeation flux improvCTuent of the helical 

module over the linear module versus Reynold's number for 2X concentration of 
whole transgenic goat milk at 298 K. 

[0038] Figure 19 is a graph of predicted permeation flux versus actual flux 

for different temporatures, wall shear rates, and milk concentrations, r2 = 0.92. 
For the helical module, the wall shear rate is given by yheUcal = 1 .95ylinear 
(Gesan-Guiziou et al., "Critical Stability Conditions in Skimmed Milk Crossflow 
Microfiltration: Impact on Operating Modes," i^/Y 80:129-140 (2000)) where 
ylinear is the wall shear rate for a linear module at the same Reynolds number, 
intemal diameter and length as the helical module. The various cases are ▲ 
Linear TS = 10%, ♦ Linear TS = 12%, ■ Linear TS = 13.2%, • Linear TS = 16.5 
to 20%, A HeUcal TS = 10%, 0 HeUcal TS = 12%, o HeUcal TS = 13.2%, n 
Helical TS = 16,5 to 20%. Wall shear rates vary from 18,500 to 87,000 s-1 and 
temperatures from 293 to 308 K. TS = total soUds. O^^otes: 1) 100% sieving is 
assumedforlactose, minerals and nonIgG whey proteins; 2) Both predicted and 
experimental flux values are based on the pressure independent values; 3) 
Electrostatic, inter-particle and particle-membrane interactions have not been 
considered.) 



DF.TATT.KD DESCRIPTION OF THE INVENTION 

[0039] The present invention relates to a method for predicting pressure 

independent praneation flux and target molecule yield in a permeate resulting 
firom crossflow membrane filtration of particles in a poly-disperse suspension. 
This method involves determining the particle size distribution of the poly- 
disperse suspension, the equivalent spherical radii of the particles, the viscosity of 
the suspension, and the maximum back-transport velocity (wO for all particles. It 
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also involves estimating the maximum aggregate packing volume fraction {^id) for 
all particles at a wall of the filtration membrane from geometric considerations, 
and selecting the particle that gives a mininium permeation flxxx: at a given 
filtration membrane shear rate, where llie selected particle has a radius (ai), and 
determining a predicted permeation flux. The method also involves determining 
packing density (^) at the membrane wall for each particle size (a, foxj ^ i) at 
the predicted peraieation flux. Also determined are the interstitial packing density 
i^mintersdce) of particlcs in the suspension which are smallest, and minimum pore 
diameter (2rnrimmum) based on the packing density of each particle. The yield of a 
target species in the filtration permeate is then estimated by calculating observed 
sieving coefficient (^o) for the target species. As a result, peraieation flux and 
target molecule yield of the poly-disperse suspension during crossflow filtration 
are predicted. 

[0040] The viscosity of the poly-disperse suspension is determined by 

experimental data or by using a modified Einstein-Smoluchowski equation: tj/ tjq = 
1+ 2.5 jfe + ki^^ where 7 is bulk fluid viscosity (kg/m.s) of the suspension, tjq is 
bulk fluid viscosity of tiie suspeosion without solute (kg/m.s), ki is particle shape 
factor (-), and ^ is particle volume fraction in the bulk suspension (-). 
Alternatively, viscosity can be determined experimentally. 
[0041] Determination of the maximum back-transport velocity («i ) is 

carried out by calculating Brownian diffusion (Jb) for all particles, where Jb = 
0A14irf^T^/?fa^Lf^^]n((^/^); calculating inertial lift (J{) for aU particles, where 
Ji = 0.036pa^//7j; and calculating shear induced diffusion (Js) for all particles 
where Js = 0.07S(a^/Lf^^jhi(^/ipO (see Table 1). ;^is wall shear rate (s"^), Tis 
temperature (K), 7 is bulk fluid viscosity (kg/m.s), ai is radius of species i(m), L is 
tube length (m), ^ is particle volume fraction at the filtration membrane (-), ^ is 
the particle volume fraction in the bulk suspension (-); and p is particle density 
(kg/m^). The maximum value ofy(Jmax) is selected for each particle. This 
maximum value of J is = t/i , the maximum back-transport for each particle in the 
suspension. 

[0042] Estimating maximum aggregate packing volxmie firaction for the 

poly-disperse suspension ((^m) at the membrane wall is carried out by determining 
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the particle sizes (ai) in fhe suspension, where at is radius of species i(m), 
detCTmining the size ratios of the particles, and using = + 0.74 (1- for a 
suspension or solution having two particle sizes such that ai>\0a2^ where is the 
maximum packing volume fraction for monodisperse spheres. Although a poly- 
disperse suspension is more complex than a susp^ision having monodisperse 
^heres, the value for a monodisperse sphere, i.e., = 0.64, can be used in the 
first iteration of the present invention in order to estimate the maximum packing 
volume fraction. Therefore, in one aspect of the present invention, ^ ~ 0.64 for 
all particles for the first iteration. In another a^ect of the present invention, 
maximum aggregate packing volume fraction at the membrane waU involves 
det^miiiing the particle sizes {a^ of species (i) in the suspension, and determining 
if the size ratio of the particles is > 10, such that > lOa,- for all a,-. If az+i > 
lOai', the maximum aggregate packing volimie fraction(^j^,;^ is calculated by 

4f« = ^ + ^ (1- 4>Mn-i\ where (pMi = <An is the maximum packing 
volume fraction for monodisperse spheres set to 0.64. 

[0043] In another aspect of the present invention, the maximum aggregate 

packing volume fraction (j^td) is estimated by determining the particle sizes in the 
suspension (a;) and calculating a maximum radius ratio of all the particles. If the 
maximum radius ratio is < 10, ^at = 0.68 is set as the value of the maximum 
packing volume fraction. 

[0044] Yet another aspect of the present invention provides that if the 

radius ratio of a suspension having two particle types is > 10, then = + 0.74 
(1- jj^), where tpm may be set to 0.64 to denote the highest packing volxmie fraction 
for a single particle. The packing volume fraction(<^ji^ at the membrane wall for a 
suspension having three particle sizes such that ai > lOaz > 100 as, is carried out 
by 

^M^^^^a- ^m) + 0.74[1 - {^n. + (1- ^)}] 

where is the maximum packing volume firaction for monodisperse spheres set 
to 0.64. An ext^sion of the above equation can be used for more than three 
particles. For a suspension where there are greater than three particle types, each 
having a size ratio greater than 10 such that, ai > lOaz > lOOas > 1000a4, etc., the 
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maximum aggregate packing volume &aciion(^Mn) at the membrane wall can be 
calculated by the recursive formula 

where ^ = is the maximum packing volume fraction for monodisperse 
spheres set to 0.64. 

[0045] Determining a predicted minimum permeation flux for a poly- 

di^erse suspension is carried out by comparing the values of (/max) determined 
for each particle type in the suspension by calculating Jb, Ju and Js as described 
above. The predicted minimum permeation flux (Jum) is determined by selecting 
from among the (/max) values for all particles the (/max) for a given particle (ai) 
having the lowest numerical value. 

[0046] The packing density {^J) for all other particles at the predicted 

permeation flux (aj for j ^ i) is also determined. This determination is carried out 
by calculating the value of (pwj such that (l>wj gives the predicted permeation flux (/) 
of selected particle (a,), xxsing the equation for back-transport that establishes 
maximum back transport for each particle (o/ for 7=/). As described above, / is 
calculated by either /b = 0. 1 lAiyj^T^lrf-c^Lf^lniif^^), /i= Qm6pa^:^li], or ^ 
/s = 0.01%{aVL)^y]xi(^^b). y is wall shear rate (s"^), Tis temperature (K), ij is 
bulk fluid viscosity (kg/!m.s), a; is radius of species z(m), L is tube length (m), ^ is 
particle volume fraction at the membrane wall (-), is the particle volume 
fraction in the bulk suspension and p is particle density (kg/m^). When / for 
(aj) has been calculated using either Jb or /s, the value of can be 'l^acked-ouf 
of the equation mathematically. However, when Ji is the governing back-transport 
mechanism, the value of caonot be directly back-calculated. In that 
embodiment of the present invention, the value of is instead calculated by 
estimating mass balance and the propensity of a particle to lift off from the 
membrane wall. In one aspect of the present invention, this can be carried out by 
determining if Uji > 10/ for the particle, where if wji > 1(X7, the determination of the 
packing density (^y) for that particle is carried out using = 0, as the particle 
back-transport velocity is much higher than the polydisperse permeation flux and 
will be readily lifted off from the membrane wall. Otherwise, i.e., if Uji < 10/, 
for that particle is carried out using 
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where j 9tjL If there is more than one particle type whose back-transport is 
governed by inertial lift, jll and jI2, their contributions to the cake at the wall can 
be approximately apportioned in the direct ratio of their volume fractions in the 
bulk suspension and inverse ratio of their back-transport by 

<*Mdn = $fejn «ji2 : (fepa Wjn 
where j ^11 or jI2 and i/jn, u^< 10/. 

This logic can be extended for more than two particle types whose back transport 
is govemed by inertial lift. 

[0047] In order to determine the permeation flux and target yield of a 

poly-disperse sxispension during crossflow microfiltration, it is also necessary to 
determine the interstitial packing density i^nnterstice) of the smallest particle in the 
suspension. This is carried out by ^wuntersuce = ^omaedliX - 2 (^jcorrected)^ where 
((Corrected = [ (<«ivi)/2 and whcTC ^ is the particle volume fraction at the 
membrane wall (-) for particle i. 

[0048] Another factor to be determined in this aspect of the present 

invention is the TniniTnum pore diameter (2rmiiiimuni) available at the membrane 
wall, based on the corrected packing density of each particle. The minimum pore 
diameter (2rnimimuni) is estimated from geometric considerations, based on a face 
centered cubic packing for the cake where th^e are four spherical particles per 
cube. Here it is assumed that the gap available for transmission in a face centered 
cubic cake is (2rminimum) = aV2 - 2ai, where a is the side of the cube and Ui is the 
radius of the particles located at the vertices and face centers of the cube. 
Therefore, in the present invention, 2rxnmimum= fli{V2[4(4/3)7iv^(*wiinterstice]^^ -2}, 
where a is radius of species / (m) and rminmmmis a minimum equivalent cake void 
radius for all cake types (m). 

[0049] The method of the present invention also involves estimating yield 

of a target molecule. This involves calculating the observed sieving coefficient 
(^q), where So = SJdl- Sa)oxp(rJ/k) + Sa). The actual sieving coefficient Sa is 
obtained from & = (*Socexp(Pem))/(5oc+exp(Pem)-l), the wall Peclet number, 
Peja is obtained from Pem = iJSJD)(Scc /s^K^X 4i is taken as the side of the face 
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centered cube of the particles of radius ai that forms the controlling cake for 
transnoission, <5bi= a = ai [(4(4/3)7t)/ ^lintersticef^y intrinsic sieving coefficient Soz is 
obtained from Soc = (1-A,)^[2 - (1- Xf] exp(-0.7146X^), X - rjrrnm. X is a ratio of 
solute to pore radii (rs/rniin)(-)> 'mm is a minimum equivalent cake void radius for 
all cake types (m), n is solute radius (m); and ^K^^ (1- X,)^'^, where is the 
equilibrium partition coefficient between membrane pore and suspension (-), and 
jfi^ is a hindrance factor for difiRisive transport (-). 

[0050] In one aspect of the present invention, the crossflow filtration 

process includes diafiltration. Diafiltration is a process whereby a filtration 
membrane is used to remove, replace, or lower the concentration of salts or 
solvents from a suspension containing biological material (Schwartz L., 
*T)iafiltration: A Fast, Efficient Method for Desalting or Buffer Exchange of 
Biological Samples," Scientific and Technical Report^ Pall Life Sciences (2003), 
which is hereby incorporated by reference in its entirety). According to the 
present invention, yield of the target species in a diafiltration experiment is 
estimated after JVa diavolumes as follows: yield =1- exp(- iVdSoavaa^)- »Soaverage is 
average observed sieving coefficient during diafiltration (-). So = 5'a)exp(- 
J/3k) + 5'a). Actual sieving coefficient iSfl is obtained from iSa = (iS'ocexp(Peni ))/(jS'oc 
+ exp(Pem ) -1), wall Peclet mmber, Pem, is obtained from Pe^ = iJSaJD)(Soc Is^ 

\ is taken as the side of the face centered cube of the particles of radius ai 
that forms the controlling cake for transmission, where ^ = a = ai 
\(^{4l2i)TC)l^iintersnc^^\ intrinsic sieving coefficient is obtained from &c = (1- 
X)^[2 - (1- Xf\ exp(-0.7146A,^), where X = rs/rnun. <^K^^ (1- ^)^'^, /-mm is minimum 
equivalent cake void radius for all cake ^es (m), rg is solute radius (m), and X is 
the ratio of solute to pore radii (ra/rnmOC-)- 

[0051] By carrying out the steps described above the pressure independent 

permeation flux and the yield of the target species are thus determined for a 
permeate resulting from crossflow membrane filtration of particles of any poly- 
disperse suspension. 

[0052] Once the predicted permeation flux and target yield are determined 

according to the present invention as described above, the performance of the 
filtration system as to permeation flux and yield can be refined using the 
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infonnation acquired by carrying out the above detenninatioiis. Thus, the present 
invention further involves re-calculating the packing density at the membrane wall 
deteraiinations for all particles in the suspension and determining if the packmg 
constraints are met for all particles. If packing constraints are not met, the 
estimations made earlier require some correction. Packing density of a particle is 
corrected by using ^icorrected = [ <^d- Th© predicted permeation flux, 

J, is then reevaluated for the particle selected as having the minimum permeation 
flux based on ^con^aed^ ^ [ The maximum back-transport velocity 

(«i), determined as described above, is also reevaluated, and ttie maximum 
aggregate packing volume fraction for all particles ((iW) at the membrane wall is 
re-estimated. 

[0053] The present invention also involves refining the determination of 

the yield of the target species. Refining the yield involves determining whether 
the suspension has a low, intracnediate, or high operating shear rate (iSo). A 
suspension is considered to have a low operating shear rate when So ^ 0.75, 
corresponding to a yield > 0.95, an intermediate operating shear rate when 0 < So > 
0.75, corresponding to yield range of from 0 to 95%, and a high operating shear 
rate when So s 0. The calculation for So is carried out using So^ SJdl- 5'a)exp(- 
J/k) + Sa) The actual sieving coefficient Sa is obtained from Sa = (S^exp(Pem 
))/(iS'oc + exp(Pem ) -1), the maximum back-transport velocity is (wOCobtained as 
described above), the wall Peclet numba:, Pcm is obtained from Pem = (JSa^/D)(S^ 
/s^Kd X <Sn is taken as the side of the face centered cube of the particles of radius 
ai that forms the controlling cake for transnoission, <Sn = ^ 

ii4(4/3)K)/^iinteiyticef^, the intrinsic sieving coefficient & is obtained from & = (1- 
Xf[2 - (1- Xf] exp(-0.7146;i^), X = rs/rmm, where A is a ratio of solute to pore radii 
(rs/rminX-), 'min is a minimum equivalent cake void radius for all cake types (m), rg 
is solute radius (m), and ^Kd^il- Xf^^ where ^ is the equilibrium partition 
coefBci©Qt between membrane pore and suspension (-), and Kd is the hindrance 
factor for diffusive transport (-). 

[0054] When an intermediate operating shear rate is determined, the value 

of (j9 is refined by calculating the stagnant film flux (J) equation for non-retentive 
membranes, where J = k In [(^ - j^enneatd) K^i - ^^enneatei)] =k]n / - iSo)], 
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wherein » <2^enneatei). So is flien corrected by replacing J= solvent 
penneation flux ( m/s) with the stagnant film flxix (J) equation for non-retentive 
membranes in the equation for observing sieving coefficient, So, where So = Sa/dl- 
5'a)exp(-J/5fc) + iSii). 

[0055] The present invention also involves constructing a plot of the 

predicted penneation flux and yield versus wall shear rate, such as shown in 
Figure lA. 

[0056] Filtration can involve noicrofiltration and ultrafiltration and may be 

carried out with a flat sheet filter, a hollow-fiber filter, or a helical filter. Suitable 
suspensions in all aspects of the present invention include, without limitation^ 
waste water, surface water, environmental pollutants, industrial waste streams, 
industrial feed streams, and streams from biomedical and bio-processuig 
industries. Such streams may contain, without limitation, proteins, cells, nucleic 
acids, colloids, milk, and suspended particles, ia any combinatioiL 
[0057] The present invention also provides a method for determining the 

packing density of particles of a poly-disperse suspension at a membrane wall. 
This method involves providing a predicted permeation flux (J)^ determining the 
packing density at a membrane wall for all particle sizes at the predicted 
penneation flux, and deternmiing the interstitial packing density {(^unterstice) of 
particles in the suspension which are smallest, thereby detennining the packing 
density at membrane wall of the particles of the poly-disperse suspension. 
[0058] Another aspect of the present invention is a method for predicting 

pressure independent permeation flux for crossflow membrane filtration of a poly- 
disperse suspension. This method involves determining viscosity of the 
suspension, detennining the maximum back-transport velocity (ui) for all 
particles, and estimating the maximum aggregate packing volume firaction {^m) for 
all particles at a wall of the filtration membrane wall from geometric 
considerations. The particle is selected that gives a minimum permeation flux at a 
given filtration membrane shear rate, where the selected particle has a radius (a/). 
A predicted permeation flux (J) is determined, and the packing density i^J) at the 
membrane wall for all particle sizes at the penneation flux (aj fory ^ i) at the 
predicted penneation flux is determined. In this aspect of the present invention 
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the determinations of viscosity of the suspension, maximum badc-transport 
velocity, and maximum aggregate packing volume firaction at the membrane wall 

for all particles are carried out using the equations given above for these 
factors. The minimum permeation flux value is selected by determining a Jmax 
value, as described above, for each particle type in the suspension, then selecting 
from all the Jmax values that Jmax having the lowest value. This method may be 
refined by also carrying out a recalculation of the packing density at the 
membrane wall determination for all particle sizes, determining if packing 
constraints are met, and correcting for packing density if packing constraints are 
not met The determinations of whether packing constraints are met or not met, 
and the calculations for correction of packing density are as described herein 
above. 

[0059] Another aspect of the present invention is a method for calculating 

yield of a target molecxile in a permeate for a poly-disperse suspension during 
crossflow membrane filtration. This method involves detennining minimum pore 
diameter (2rniinimum) based on the packing density of each particle, and estimating 
yield of a target species in the filtration permeate by calculating observed sieving 
coefficient (^o) for the target species. Determination of the minim um pore 
diameter (2rniiiiiiiiuni) is carried out as described above. The determination of 
minimum pore diameter (2rniinimuin) is then used to estimate the yield of a desired 
target molecule in the poly-disperse suspension by carrying out the calculation for 
the observed sieving coefficient (iSo), as described above. This aspect of the 
present invention may fiirfher involve diafiltration. The yield of the target species 
on a diafiltration experiment can be estimated after JVd diavolumes, as described 
above. 

[00601 Th^ present invention also relates to a method for designing a 

crossflow membrane filtration system for a poly-disperse suspension. The 
p^ormance parameters of a crossflow membrane filtration system can be 
designed for any selected poly-disperse suspension by applying the methods 
described herein for predicting pressure independent praneation flux and 
detennining target molecule yield. This involves determining the particle size 
distribution of the poly-disperse suspension and the equivalent spherical radii of 
the particles, and detennining the viscosity of the suspension and the maximum 
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back-transport velocity (^0 for all particles. It also involves estimating the 
maviTniiTTi aggregate pacldng volume fraction (^Af) for all particles at the filtration 
mCTibrane from geometric considerations; selecting the particle that gives a 
TTiiniTnnm permeation flux at a given filtration membrane shear rate, where the 
selected particle has a radius (aO, and detenniiiing a predicted permeation flux. 
The metiiod also involves determining packing density at the membrane wall for 
all particle sizes at the predicted permeation flux, interstitial packing density 
i^intenstice) of particlcs in the suspension which are smallest, and minimiTm pore 
diameter (ZrndimniiiiO based on the packing density at the membrane wall of each 
particle. The yield of a target species in the permeate is then estimated by 
calculating observed sieving coefificient (So) for the target species. All of these 
determinations are carried out as described above for other aspects of the present 
inventioiL Conditions for filtration based on the prediction of permeation flux and 
target molecule yield are then optinaized to design a filtration system for the 
selected poly-disperse suspension. 

[0061] Yet another aspect of the present invention is a method of selecting 

operating conditions of a crossflow filtration system for poly-disperse 
suspensions. This method involves predicting the pressure independent 
permeation flux and target molecule yield in a permeate resiilting from crossflow 
membrane filtration of particles in a poly-disperse suspension as described in 
detail above. This predicts permeation flux (process time) and target molecvile 
yield of the poly-disperse suspension during crossflow membrane filtration. 
Operating conditions of the system are selected using the determination of 
limiting pressure independent permeation flx3x for a given shear rate to obtain an 
optimal balance between permeation flux and yield of a target species. 
[0062] The present invention also relates to a method of modeling a 

process for filtration of a poly-disperse suspension. This method involves 
applying the method for predicting pressure independent permeation flux and 
target molecule yield in a permeate resulting from crossflow membrane filtration 
of particles in a poly-disperse suspension, using the calculations described above, 
and using a computer-generated program to model the process for filtration of a 
poly-disperse suspension. 
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EXAMPLES 
Example 1 - Aggregate. Transport Model 

[0063] With the framework described above, the aggregate trajosport 

model of the present invention is summarized below. Refer to Table 2 for a 
definition of symbols used herein* 

[0064] Step 1 . Determine the particle size distribution of the feed 

suspension and evaluate the equivalent spherical radii. This can be obtained from 
literature, by size exclusion chromatography or by membrane firactionation. 
[0065] Step 2. Evaluate the viscosity of the suspension by experiment or 

estimate it using the modified Einstein-Smoluchowski equation (1) 

77/770 = 1+ 2.5^ + ki<zJb^ (6) 

[0066] Step 3, Evalxiate the maximum back-transport velocity, «i, for a 

particle based on Brownian difEusion, shear induced difiiision, and inertial Uft at 
the proposed operating wall shear rate assxmiing full retention for all solutes, using 

Mi = Max[B(ai, y, <j>^ (^,L, 77), S(ai, /, L, 7i\l{aur. 7)] (J) 

where B, S, and I denote the fimctionaUties for Brownian diffusion, shear induced 
difiiision and inertial lift models, respectively, shown in Table 1 , above. It is 
suitable to set ^ = 0.64 for each species for the first iteration ^odds, "The 
Porosity and Contact Points in Multicomponent Random Sphere Packings 
Calculated by a Simple Statistical Geometric Model," J. Colloid Interface 77:317- 
327 (1980), which is hereby incorporated by reference in its entirety). 

[0067] Step 4. Estimate the maximxim aggregate packing volxmie fraction 

for all particles, ^ at the wall, from geometric considerations. For the poly 
disperse case, this could be much larger tiian the widely used value 0.64 
depending on the size ratios of the particles. If the size ratio is more than 1 0, the 
small particles behave as a continuous fluid with respect to the large particles and 
can migrate into the interstices easily (Farris, *Trediction of the Viscosity of 
Multimodal Suspensions from Unimodal Viscosity Data," Trans. Soc. Rheol. 
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12:281-301 (1968); Piobstein et al., '"Bimodal Model of Concentrated Su^ension 
Viscosity for Distributed Particle Sizes," J, Rheol 38:81 1- 829 (1994); Gondret et 
al., *T)ynanuc Viscosity of Macroscopic Suspeasions of Bimodal Sized Solid 
Spheres," J. Rheol 41:1261-1274 (1997), which are hereby incorporated by 
reference in their entirety). For example, for a poly-disperse mixture comprising 
particles of three sizes such that a\ >\Qa2 > lOOas the relation (8) may be used 

^W=««b+jAn(l-iU + 0.74[1 -{^+(f^ (1- $^b)}] (8) 

where ^ is the maximum packing volume fraction for monodisperse spheres = 
0.64 (36). In this special case, ^ = 0.96. 

[0068] Step 5. Repeat step 1 for all particle sizes and select the particle 

that gives the miiiimum permeation flux at the given wall shear rate. The 
corresponding permeation flux is the predicted one, 

J= Min [w, M2,. . -t^n]. The selected particle has a radius ai (9) 

[0069] Step 6. Evaluate packing density for other particle sizes (aj for j ^ 

i) at this permeation flux. Calculate ^ from the equation 

J = Max[ B(aj, A 7). S(aj, ^, L, rj), I(aj, ^ ^, X, tj) ] (10) 

for all j using the equations for J shown in Table 1 . For particles whose back 
transport is governed by inertial lifl; [Table 1, equation (3)], the wall concentration 
cannot be calculated explicitiy by equation (10), but can be estimated by using 
mass balance and the propensity to lift off. For example, if there is only one such 
particle type, jl, check if mji > 10 J. If yes, set = 0 as the particle back-transport 
velocity is much higher than the polydisperse permeation flux and will be readily 
lifted off from the membrane wall. If Uji is <10J where one particle type is 
govemed by inertial lift, then is calculated by 

^1= ^M'^ (11) 
whCTe j 9£jL If there is more than one particle type whose back-transport is 
govemed by inertial lift, jll and jI2, their contributions to the cake at the wall can 
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be approximately apportioned in the direct ratio of their volume fractions in the 
bulk solution and inverse ratio of their back-transport by 

J^vsrjii + <hin = i^ji/ - S (12) 
<*wjn • = ^ii2 • (^ji2 Wjii (13) 
where j ^gll or jI2 and uju, Wji2< 10/. 

This logic can be extended for more than two particle types whose back transport 
is governed by inertial lift 

[0070] Step 7. Check 2 fJwi < and other packing constraints. These 

depend on the particle sizes in the cake and have to be developed specifically for 
each case. In Example 1, below, a typical case of a tridisperse suspension with 
two large and one small particle type is illustrated. The packing constraints for 
this case are depicted in Figure IC, Guidelines for developing packing constraints 
are described in Example 3, below. If packing constraints are satisfied, go to Step 
8, or else use 

Corrected ^ ^ [ (^)/2 ^ (14) 

For the particle selected in Step 5, re-evaluate J based on ^wicorreaed instead of 0.64 
by repeating Steps 3 and 5. 

[0071] Steps. Evaluate interstitial packing density, ^//wteryftce of the 

smallest particle by 

^iinterstice ~ ^mcorrectedl 

(1 - S ^corrected) (15) 

[0072] Step 9. Based on the corrected packing density of each particle, 

estimate the minitmiTn pore diameter 2rnimimum from geometric considerations 

2r„unimum = ai{V2[4(4/3)7i/^J^^>^^„^ J 1/3 -2} (16) 

This is based on a face centered cubic packing for the cake where there are four 
spherical particles per cube. Here it is assumed that the gap available for 
transmission in a face centered cubic cake is 2rminimum = «'>/2 — 1a\ where a is the 
side of the cube and Ui is the radius of the particles located at the vertices and face 
centers of the cube, a is evaluated based on the value of and ^onterstice- 

If this 

gap is greater than 2rs, a particle of radius rs can be transmitted through the cube 
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by displacing two particles located at face centers on the way. Please refer to 
Example 2 and Figure 1 C for further details. 

[00731 Step 10, Estimate the yield of target species (with equivalent 

radius = ts) in the permeate by calculatiag the observed sieving coefficient. So 
(30). 

So = SaKil- Sa)CXpi-J/k) + Sa ) (17) 

where actual sieving coefficient Sa is obtained from (30) 

Sa = (Soc&qp(Peu, )y(S^ + exp(Pein ) -1). (18) 

The wall Peclet number, Pcjn is obtained from 

Pen, = (JSJD)(So. /s^Ka ) (19) 

where is taken as the side of the face centered cube of the particles of radius ai 
that forms the controlling cake for transmission. In general, the governing case 
for flux and product transmission (corresponding to rminimum ) may be different. 
Hence, 

• Sn, = a = ai[i4i4/3)K)/i^iintersacef^ (20) 

The intrinsic sieving coefficient »Soc is obtained from (40) 

So, = (l-;i)^[2 - (1- Jlf] exp(-0.7146;i^) (21) 
where il^rs/rmimimiiii. is estimated from (41) 

^Ka = il^Jlf^ (22) 

[0074] The yield of the target species in a diafiltration experiment can be 

evaluated after JVd diavolumes by using the following relation (30) 

Yield = 1- exp(- iVk^oaveragc) (23) 

[0075] Step 1 1 . There are three possible scenarios corresponding to low, 

intermediate and high operating shear rates. For low operating shear rates, the 
cake will be dominated by the larger particles leading to high observed sieving 
coefficients where So ^ 0.75, corresponding to Yield > 0.95 according to equation 
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(20) for 4 diavolumes. For such cases no further refinement is needed. For 
intermediate operating shear rates, 0<So< 0.75, leading to a yield range from 0 to 
95%. For this case, is further corrected by using the stagnant film flux equation 
for non-retentive membranes (for ^ » ^pame&td) 

7= hi - ^ermeatd) K^i ' ^^nneatei)] ^k]n[^ /(fei(l - So)] (24) 

for the transmitted species in Step 3. Steps 3 to 10 are repeated until the values of 
So obtained by equations (14) and (20) are within 10% of each other. For high 
operating shear rates, the cake is dominated by the smallest particles and So = 0 
implying no transmission of the smallest particles. Thus, in tibis case, the original 
assumption of full retention is valid and no further iterations are required. 
[0076] Step 12. Construct plot of predicted permeation flux and yield 

versus wall shear rate for the pressure-independent regime. 
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Table! 



Symbol 


Refers To: 


a 


Radius of species (m) 


D 


Molecular diffusion coefficient (m^/s) 


J 


Solvent permeation flux ( m/s) 


k 


Mass transfer coefficient (m/s) 


h 


Shape factor (-) 




Hindrance factor for diffusive transport (-) 


L 


Tube lengfh, (m) 


In 


Natural log 


iVd 


Niimber of diavolumes during diafiltration (-) 


Pem 


Membrane Peclet number (-) 


n 


Equivalent cake void radius (m) 


^ min 


Minimum equivalent cake void radius for all cake types (m) 




Solute radius (m) 


So 


Observed sieving coefficient (-) 


'Soaveraee 


Average observed sieving coefifidlent during diafiltration (-) 


Sa 


Actual sieving coefficient (-) 




Asymptotic (intrinsic) sieving coefficient (-) 


T 


Temperature (K) 


Ui 


Back dif&sion velocity of particle i (m/s) 


s 


Momentum boundary layer thickness (m) 




Membrane/cake thickness (m) 


6 


Cake/membrane porosity (-) 


^. 


Equilibrium partition coefficient between membrane pore 




Particle volume firaction in the bulk suspension (-) 




Particle volume fraction at the membrane wall (-) 


Y 


Wall shear rate (s'^) 


V' 


Bulk fluid viscosity (kg/m.s) 




Bulk fluid viscosity without solute (kg/m.s) 


K 


Boltzmaim constant (J/mol K) (1.3807 x 10"^ J K"') 




Ratio of solute to pore radii (rs/rmin)(-) 


P- 


Particle density (kg/m^) 
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F vam ple 2 - Sample Calculation of Permeation Flux and Yield for a 
Hypotiietical Poly-disperse Suspension 

[0077] The method of the present invention is first illustrated by applying 

it to a hypothetical suspension comprising three diffarraat sized particles: 1 0, 1 80, 
and 300 nm. The relevant data for the tridispCTse suspension and microfiltration 
system are shown below in Table 3 and the calculations are descaibed step-wise 
below. 

TABLES 



Module length - 


0.3 m 


Number of fibers 


6 


Fiber LD. 


1,27 nmi 


Pore size 


0.1 \im 


Temperature 


298 K 


Shear rate 


32,400 sec'" 


Particle sizes and bulk 
concentrations 


lOmnat 1%, 
180 nm at 5%, 
300 nm at 6% by 
volxime 


Target particle radius 


lOnm 



[0078] Calculations were carried out as follows. 

[0079J Step 1 : Evaluate particle size and equivalent radii. Particle sizes 

are given above in Table 3. 

[0080] Steps 2-7: Evaluate the pressure independent permeation flux at the 

given conditions. This is evaluated iteratively by a stpreadsheet program 
(incorporating Steps 2-7 of the method of the present invention) at the sheax rate 
of 32,400 sec'*. J= 30 Imh and the selected particle has a radius = 180 nm and 
A-isonm = 0.435. In this example, there are three different particle radii, 10 nm, 
1 80 run, and 300 mn. The maximum packing firaction considering the 1 80 nm and 
300 nm will be 0.68 (Gondret et al., 'T)ynamic Viscosity of Macroscopic 
Suspensions of Bimodal Sized SoUd Spheres," J. Rheol. 41:1261-1274 (1997), 
which is hereby incorporated by reference in its entirety). The 10 nm particles are 
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less than 1/lOth of the next higher size particles and hence, can freely move 
aroxmd the interstices of the bigger particles (Probstein et al., '"Bimodal Model of 
Concentrated Suspension Viscosity for Distributed Particle Sizes," J. Rheol 
38:811-829 (1994), which is hereby incorporated by reference in its entirety). 
Assuming that the 1 0 nm particles can reach a maximum pacldng of 0.74 
(corresponding to face centered ciibic packing) in the interstices of the bigger 
particles, the maximum aggregate packing = 0.68 + 0.74(1-0.68) = 0.92. This 
value is used for ^ia steps 4 and 7. 

[0081] Step 6: Evaluate the packing densities of the other particles. The 

spreadsheet program based on Step 6 gives fJwionm = 0,27 and jJwaoonm = 0.179. 
[0082] Step 8: Evaluate the interstitial packing density of the smallest 

particle. ^lOnm interstice ^ 

0.27/[l-(0.435 + 0.179)] = 0.7 using equation (12). 
[0083] Step 9: Evaluate the minimum pore radius for transmission. Using 

equation (13), 10.38 nm. 

[0084] Step 10: Estimate the yield of the target species. Using % = 

rs/rndnimum, A, = 0.963, Usiug equations (18) and (19), ^'oc = 0.001382 and ^Kd = 
3.446E -7, where Pem = (•/4x/D)(&c /s^K^ )(16). Cake thickness, = 28.8 nm 
using equation (20), s = 1- ^lotsm interstice = 0.3. Next, D is estimated. Using 
Einstein's equation, D = (tdT)/(6n7] at) = 1.67x 10'^^ m^/sec. Using the value of J 
calculated earlier, P^m = 384. Using equation (16) and (17), 5'^= 0.00138 and So = 
0.088. Finally, using equation (20) for iSoaverage - 0.088 and for a diafiltration run 
after 4 diavolumes, yield is determiaed to be 30%. Therefore, for this sample 
problem, J— 30 Imh and Yield = 30% for the first iteration. 
[0085] Step 1 1 : RejSne the yield of the target species. Adopting 

*hl[(<j)wi-4>permeateO/(<l>W-(|>p™^ (24) 

for the 10 nm particles as 0 < 5i> < 0.75 Step 9 is executed to obtain, after a few 
iterations, 31,5 hnh, i^tvionm = 0.205, j^isomn = 0.49 and ^Jwsoonm = 0.19. 
Following the same procedure as for the jSrst iteration. So = 0.37 and after 4 
diavolumes for 10 ran particles, assuming iSoaverage = 0.37 the Yield = 77%. 
Therefore, for this sample problem, J= 3 1 .5 Imh and Yield = 77%. 
[0086] On the other hand, if the shear rate was increased to 44,600 sec"^, 

the pressure-independent permeation flux would have increased to 36 hnh but the 



i 
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yield of the 10 mn radius target species is predicted to be 0! In actual practice, 
there will be a small finite yield because close to the exit of fbe module the 
permeation flux will be less than the pressure-independent permeation flux due to 
low transmembrane pressure in that region. The model is expected to be more 
accurate for the uniform transmembrane pressure case or in the case of a stirred 
cell where the transmembrane pressure is constant. 

[0087] The sample calculations for a hypothetical suspension comprising 

particles of three sizes indicate three dijSor^t curves for pressure-independent 
permeation fltix versus shear rate. This is shown in Figure 2A. As shown in 
Figure 2A, at high shear rates the smallest particles have the lowest back-transport 
velocities. Hence, at high shear rates the poly-disperse flux will be determined by 
the smallest particles. 

[0088] The predicted poly-disperse permeation flux and packing densities 

in the cake for each particle size at different shear rates are derived and plotted as 
shown in Figure 2B. It is observed that the packing densities of large 300 and 180 
nm particles decrease with increasing shear rate making way for the smallest 
particle size of 10 nm. The transition from a coarse to fine cake occurs at a shear 
rate of around 40,000 s'\ Operation at shear rates lower than 40,000 s"^ will lead 
to higher transmission of 10 nm particles because of the coarse nature of the cake. 
The converse is expected at shear rates higher than 40,000 s"\ This methodology 
ignores inter-particle interactions. At low shear rates the packing d^isities of the 
larger particles tend to increase. When applying this model to real situations it has 
to be determined if the preponderance of larger particles in the cake lead to 
deposits of smaller particles due to adhesion, hydrophobic interaction, etc. If that 
is tihe case, the best solution may be to operate at high shear rates but, in the 
pressure dependent regime, which will lead to a sparse cake of fine particles 
allowing good transmission of the target particles. 

[0089] The sensitivity of permeation flux with respect to module length 

and total sohds volume fraction in the bulk suspension is shown in Figure 3 and 
Figure 4, respectively. As expected, the model predicts higher permeation fluxes 
with shorter module lengths and for feeds of lower bulk concentrations. Figure 5 
depicts the expected yield (defined as the fraction of target molecules originally in 
the feed that is recov^^ in the permeate) of the 10 nm radius particles against the 
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wall shear rate in the pressure indepCTid^t permeation flux regime. The yield is 
close to 100% until a shear rate of - 28,000 s'^ is reached. Beyond this value of 
tiie shear rate, the yield drops rapidly. This captures the phenomenon reported by 
sev^ul researchers (Bene et al., "Skim Milk Crossflow Microfiltration 
Performance Versus Permeation Fl\ix to Wall Shear Stress Ratio," J. Membr. ScL 
1 17:261-270 (1996); Gesan-Guiziou et al., "Critical StabiUty Conditions in 
Crossflow Microfiltration of Skimmed Milk: Transition to Irreversible 
Deposition,' V. Membr, ScL 158:211-222 (1999); Gesan-Guiziou et al., 
*Terfonnance of Whey Crossflow Microfiltration During Transient and Stationary 
Operating Conditions," Membr. Set 104:271-281 (1995); Gesan-Guiziou et al,, 
"Critical Stability Conditions in Skimmed Milk Crossflow Microfiltration: Lnpact 
on Operating Modes," Lait 80:129-140 (2000), which are hereby incorporated by 
reference in their entirety) of a critical permeation flux to shear rate ratio beyond 
which protein transmission during microfiltration of milk reduces drastically. 
Thus, this model can be used to specify the limiting pressure independent 
permeation flux for a given shear rate and the expected yield of a target protein 
under these conditions. To obtain good yield, one should opiate at much lower 
permeation fluxes. Depending on the size of tiie target species, an optimization 
exercise can be carried out in order to choose a desirable balance between the 
yield of the target species and the permeation flux, hi typical hollow fiber 
modules, the transmembrane pressure cannot be varied independentiy of the shear 
rate because a high shear rate gives rise to a high axial pressure drop and, 
consequentiy, a high transmembrane pressure, high permeation flux, and, 
possibly, a lower yield. This leaves the option of reducing p^meation flux as 
required, by throttling the permeate stream. This, however, leads to inefficient 
usage of the module as there will be reverse permeation called Starling flow 
(Hammer et al., "Quantitative Flow Measurements in Bioreactors by Nuclear 
Magnetic Resonance Imaging," Bio/Technology 8:327-330 (1990); Heath et al., 
^•Magnetic Resonance hnagmg and Modeling of Flow in Hollow-Fiber 
Bioxeactor,''' JIChE Journal 36(4):547-558 (1990), which are hereby incorporated 
by reference in their entirety) at downstream locations in the module where the 
permeate pressure will exceed the tube pressure. This is depicted in Figure 6 A. 
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[0090] In addition, even at low transmembiane pressures, a compact cake 

may form near the entrance of the filter due to higher transmembrane pressure and 
resultant penneation flux in this region. Therefore, good sieving of the target 
molecule will only be achieved after some distance along the membrane. To 
alleviate this difficulty, short modules to reduce Ihe axial pressure drop are 
recommended, as shown in Figure 6B. Better still, this problem is completely 
eliminated by using pameate re-circulation to maintain uniform transmembrane 
pressure along the len^ of the module, as shown in Figure 6C. In this way, hi^ 
shear rate operation with low permeation is possible and high yields of target 
species can be obtained. 

[00911 Another beneficial mode of operation can be to maintain the wall 

volume firaction ^ of the speraes corresponding to the least void radius, constant 
as evaluated in step 6 of the method. This will result in a feirly uniform sieving 
coefficient during diafiltration. To achieve this, the constant <p^, the method 
recommended by van Reis et al (Reis et al., "Constant Cwau Ultrafiltration Process 
Control," J. Membr. ScL, 130:123-140 (1997), which is hereby incorporated by 
referoice in its entirety) could be adopted. 

Example 3 - Guidelines for Developing Packii^ Constraints 

[0092] Packing constraints of the cake formed at the membrane waU 

depend on the size distribution of the particles in the bulk suspension. A few 
aspects have been covered in Steps 4 of the aggregate transport model described 
above herein, and as shown in Figure IC. Further guidelines to evolve packing 
constraints are presented here. 

[0093] StepBl. Estimate the maximum aggregate packing volume 

fiaction for all particles. Variants of equation (8) of Step 4 may be used. If the 
maximum radius ratio of the particles is < 10, can be set to 0.68 based on 
reference (Gondret et al., •T)ynamic Viscosity of Macroscopic Suspmsions of 
Bimodal Sized SoUd Spheres, J. Bheol. 41:1261-1274 (1997), which is hereby 
incorporated by reference in its entirety). If there are two distinct groups of 
particles separated by a fector of 2: 10 in radii, a truncated version of equation (8) 
may be used 
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^W=iftn + 0.74(l-^!b) (Bl) 
where (pm may be set to 0.64 to denote tiie highest packing volume fraction for a 
single species. In a manner similar to equations (8) and (Bl), ^ii/ for the case for 
more than three distinct particle size groups can be estimated. The particle 
composition of the cake and the bulk suspension will be different because of the 
different back-transport mechanisms applicable for different particle types. It is 
possible that certain particles get swept away from the wall at very high back- 
transport rates. These particles can be eliminated from the cake if their back- 
transport rates are more than 10 times higher than the poly-disperse flux evaluated 
in Step 5. This will simplify the problem greatly. 

[0094] Step B2. Evaluate the interstitial packing of the smallest particles. 

It is assiuned that the smallest particles can pack tightiy in an FCC structure 
within the interstices of particles which have radii more than 1 0 times that of the 
smallest particles. If there is a range of sizes within the smallest size group, the 
mean radius may be considered as an approximation. Thus: 



^smaUest<0^1'^{l'^</^]) (B2) 

^snuaiest^O^eA (B3) 

^smallest interstice = {^smalles^f (1" ^^j) C^) 

where j ^ smallest particle. 



Example 4 - Gap Calciilations for a Face Centered Cubic Arrangement of 
Particles 

[0095] hi a face centered cubic arrangement each face center particle is 

shared by two cubes and each comer particle is shared by eight cubes. Hence, 
total number of particles per cube = 6/2 +8/8 = 4. Let a = side of the cube, a,- = 
radius of the particle and ^untersoce = interstitial packing fraction of the particle, 
then 

i^ninte^tice = (4(4/3)7C a,- (CI) 

After rearrangement, 

a = at [(4(4/3)71)/ ^imtersticei^ (C2) 
The gap for particle transmission 

aV2 - 2ai = ai{M'^{4l3)nli^iintersticef^ -2} (C3) 
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Example 5 - Application to a Model for Microfiitratioii of IgG from 
Transgenic Goat Milk 

[0096] The mefliod of the present invention was used to design a 

predictive aggregate transport model to meet the technical challenge of recovering 
human IgG fusion protein from transgenic whole goat milk by microfiltration at 
reasonable cost with high purity and yield. To test the model's predictability of 
permeate flux and mass transport, a comprehensive series of experiments with 
varying wall shear rate, feed temperature, feed concentration, and module design 
are presented herein. 

[0097] A very good fit was obtained between the model predictions and 

measurements for a wide variety of experimental conditions. For microfiltration 
module design comparison, a linear hollow fiber module (representing current 
commercial technologies) gave lower permeation flux and higher yield than a 
helical hollow fiber module (represmting the latest self-cleaning methodology). 
These results are easily explained with the model which is now being used to 
define operating conditions for maximizing performance. 

[0098] The procedure described by the model is generally applicable and 

can be used to obtain optimal filtration performance for applications other than 
milk. 

[0099] Here, the aggregate transport model of the present invention is 

tested using whole transgenic goat milk, an enormously complex and challenging 
fluid, for recovery of a desirable molecule such as a heterologous immunoglobulin 
(IgG). The transg^c process has evolved recenfly as an economically attractive 
way of producing large amounts of human therapeutic proteins (Kreeger, 
•Transgenic Manomals Likely to Transform Drag Making,*' The Scientist 1 1(15) 
1997); Pollock et al., 'Transgenic Milk as a Method For The Production of 
Recombinant Antibodies," J. Immunol Methods 231:1.47-157 (1999); John et al., 
"E3q)ression of an Engineered Form of Recombinant Procollagen in Mouse Milk," 
Nature Biotech. 17: 385-389 (1999); which are hCTeby incorporated by reference 
in their entirety). In contrast to the traditional method of using large scale batch 
cell cultures or blood plasma fractionation, transgenic production involves flie 
creation of genetically altered animals which express the desired protein in their 
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millc. A DNA construct, comprising the sequence that will encode the target 
hximan protein and an adjacent pronoLOter sequence which facilitates expression 
only in the mammary glands, is inserted into a goat cell line by transfection. The 
nucleus is removed from an oocyte which is extracted from an animal. A 
transfected, selected transgenic cell is then ftised with the enucleated oocyte by 
electrofusion. After 24-48 hours in culture, the embryo is transferred to a 
surrogate mother. The putative transgenic animals are identified by screening the 
oflfepring for the transgene by PGR and Southern blotting. After the selected 
females mature, they are bred and the milk produced after gestation is tested for 
protein expression. The process therefore involves two gestation periods and one 
maturing period. For goats and cows this period is 16-18 months and 3 years, 
respectively OPollock et al., 'Transgenic Milk as a Method For The Production of 
Recombinant Antibodies," J. Immunol. Methods 231:147-157 (1999), which is 
hereby iucorporated by reference in its entirety). 

[00100] Several companies and organizations now use this technology. A 
plethora of therapeutic proteins such as human monoclonal antibodies, tissue 
plasminogen activator, antithrombin III, and human lactoferrin are proposed for 
manufacture by the transgenic process and are in various stages of FDA approval 
(Kreeger, "Transgenic Mammals Likely to Transform Drug Making/' Tlxe 
Scientist 1 1(15) 1997); PoUock et al., ^Transgenic Milk as a Method For The 
Production of Recombinant Antibodies," J. Immunol Methods ^231:147 -157 
(1999); Prunkard et al,, "High-Level Expression of Recombinant Hxraian 
Fibrinogen in the Milk of Transgenic Mice," JViz^Mre^io^^^ 14: 867-871 (1996); 
Mckee et al., *Troduction of Biologically Active Salmon Calcitonin in the Milk of 
Traiisgenic Rabbits," JVaft/re5/o^ec/i, 16:647-651 (1998); John etal., "Expression 
of an Engineered Form of Recombinant Procollagen in Mouse Milk," Nature 
Biotech. 17: 385-389 (1999), which are hereby incorporated by reference in their 
entirety). Remaining technical challenges include proper post translational 
modification of the secreted proteins and adequate product recovery from milk. 
[00101] Although milk from transgenic farm animals can become a large 
source of therapeutic proteins, the complexity of milk combined with the low 
concentration of target protein compUcates the recovery process. Whole milk 
consists of more than 100,000 different molecules disposed in three phases 



wo 2004/016334 



-40- 



PCTAJS2003/025230 



namely, lipid casein, and whev ( Dairy Processing Handbook. TetraPak 
Processing Systems, AB, S-221 86, Lund Sweden, (1995), which is hereby 
incorporated by reference in its entirety). The composition and properties of the 
main constituents in goat milk are given below in Table 4. Essentially, goat milk 
consists of 4 wt % each of protdn, fat, and low molecular weight moieties like 
carbohydrates, sugars, and salts. About 80 % of the proteins exist in the form of 
casein micelles. Heterologoxis recombinant proteins can be overproduced in the 
range of 0.2 to 1 wt %. The first stqp in isolating heterologous proteins fi:om 
transgenic milk involves the removal of casein micelles and fat globules firom the 
milk leaving behind low molecular weigiht salts and sugars. The traditional 
methods used by the dairy industry to isolate proteins from milk include 
pasteurization followed by enzymatic coagulation or add precipitation at pH 4.6 
(pi of casein). These steps are often unsuitable for the recovery of heterologous 
proteins, because they can be temperature and pH sensitive. Additionally, the 
coagulation process traps most of the target protein within casein pellets resulting 
in poor yields (Morcol et al., **Model Process for Removal of Caseins from Milk 
of Transgenic Animals," jB/otecAno/. Prog. 17:577-582 (2001), which is hereby 
incorporated by reference in its entirety). The removal of fat is another important 
issue. There are inherent difficulties with centrifugation regarding scale-up and 
contamiaation. Milk is skimmed industrially at a centrifugal force of 500g. Fat 
globules smaller than 5 jmi escape this process. Sub-micron size fat globules can 
be removed by ultracentrifugation (600,00g), but this is only practical for small 
volumes (Gardner et al., "Delipidation Treatments for Large Scale Protein 
Purification Processing," Master's thesis at Virginia Polytechnic Institute and 
State University (1998), which is hereby incorporated by referaice in its entirety). 
Instead, microfiltration can be used for the removal of casein and fat (retained in 
the retentate) with the target protein passing with the permeate (Pollock et al., 
"Transgenic Milk as a Method For The Production of Recombinant Antibodies," 
y. Immunol Methods 231:147-157 (1999); Meade et al.. Gene Expression 
Systems: Using Nature for the Art of Expression, Academic Press, pp. 399-427 
(1 999); which are hereby incorporated by reference in their entirety). Thus, 
microfiltration followed by ultrafiltration with various chromatographic steps 
becomes an attractive method for transgenic noilk processing. Many of these 
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processes have been patented for similar applications (U.S. Patent No. 5,756,687 
to Denman et al.; U.S. Patent No. 6,183,803 to Morcol et al., which are hereby 
incorporated by reference in their entirety). Transgenic milk is neither pasteurized 
nor homogenized in order to prevent damage and loss of the target heterologous 
proteins. Fat globules and casein micelles are the putative foulants for whole milk 
microfiltration, because the other moieties are much smaller than the average pore 
size of the 0. 1 jim microfiltration membrane, they easily pass through the 
membrane with the permeate. 

[00102] Most of the lipids (> 95%) in milk exist in the form of globules 
ranging from 0.1 to 20 ^tm in diameter (Gofif et al., "Dairy CSaemistry and 
Physics," in Dairy Science and Technology Handbook, Vol. 1, Principles and 
Properties, Hui, Y.H. ed.. New YorkiVCH Publishers, Chap. 1, pp. 1-81 (1993), 
which is hereby incorporated by reference in its entirety). In non-homogenized 
milk, the liquid fat droplets are encased by a 8 to 10 nm thick membrane called 
the native fat globule membrane (FGM). This is comprised of apical plasma 
membrane of the secretory cell which continually envelopes the fat droplets as 
they pass into the Ivimen of the mammary gland. The FGM is therefore, composed 
of phospholipids and proteins and is characterized by a very low interfacial 
surface tension, 1 to 2.5 mN/m, between the fat globules and the serum phase. 
This prevents the globules from floccxdating and from enzymatic degradation. 
Homogenization decreases the diameter of the fat globules, thereby significantiy 
increasing the surface area of the fat globules resulting in insufficient native FGM 
to cover all the fat globules. On disn^tion of the native FGM, the interfacial 
surface tension increases to a value of '--15 niN/m allowing sermn proteins and 
casein micelles fireely to adsorb onto the exposed fat globules. Thus, 
homogenization leads to a reduction in the size of fat globules as well as loss of 
proteins through adsorption (Meade et al., **Gene Expression Systems: Using 
Nature for the Art of Expression," Academic Press, pp. 399-427 (1999), which is 
hereby incorporated by refereaice in its entirety). The latter effect is expected to 
reduce the yield of target protein and the former effect increases membrane 
fouling because of a lower value of back-transport due to shear or inertial lift 
Thus, the large fat globules in non-homogenized whole milk may not be the chief 
foulant in whole milk microfiltration. This hypothesis is tested herein. 
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[00103] There is considerable dispute as to the exact nature of the casein 
micelle. The present consensus is that tiie casein micelle is a roughly spherical, 
fairly swollen particle of 0.1 to 0.3 pm diameter with a hairy outer layer (Walstra, 
"Casein Sub-Micelle: Do They Exist? Int Dairy J. 9:189-192 (1999), which is 
hereby incorporated by reference in its entirety). This is supported by electron 
microscopy studies (McMahon et al., "Rethinking Casein Micelle Structure Using 
Electron Microscopy," J. Dairy Set 8 1 :2985-2993 (1 998), which is hereby 
incorporated by refermce in its entirety). The hairy layer is comprised of C- 
tenninal ends of K-casein. TMs prevents furflier aggregation of imcelles and 
flocculation by steric and electrostatic repulsion at pH values higher than 4.6, the 
pi of casein. Thus, at the physiological pH of milk, i.e., 6.4-6.6, the casein 
nacelles predominantly exist as distinct particles of a size range comparable to the 
mean pore size (0. 1 ^ixn) of the poly(ether sulfone) microfiltration membrane. 
This is expected to result in a low sbear-induced diffusion coefficient as well as 
fouling by pore blockage, cake formation, and pore constriction for larger pores. 
It is thus expected that for this case, the casein micelles are the main candidates 
for pore plugging and cake formation (fouling). This is corroborated by 
polyacrylamide gel electrophoresis studies of permeate samples of m il k clarified 
by microfiltration with a 0.2 |jm average pore size ceramic membrane which 
indicate negligible casein transmission through the membrane. 
[00104] Here, crossflow microfiltration of raw goat milk is carried out, the 

first step in the protein recovery process firom transgenic whole goat milk, A 
working predictive model for describing the rather complex process of transgenic 
whole milk microfiltration has been developed applying the method of the present 
invention. The next step is to develop an optimizing strategy for diafiltration 
using this model and theai to conduct diafiltration experiments as a validation. In 
this study experiments were undertaken for both linear and helical hollow fiber 
modules (Luque et al.^ "A New Coiled Hollow-Fiber Module Design for 
Enhanced Microfiltration Performance in Biotechnology," Biotechnol Bioeng. 
65:247-257 (1999), which is hereby incorporated by reference in its entirety). It is 
well known that flie main problems in milk microfiltration are low flux and poor 
protein transmission due to concentration polarization and fouling. High 
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cTOSsflow velocity, back pulsing, ptilsatile flow, and Taylor and Dean vortices 
represent some of the techniques used by researchers to mitigate this problem 
(Belfort et al., •The Behavior of Suspensions and Macromolecular Solutions in 
Crossflow Microfiltration," J. Membr. Sci. 96:1-58 (1994), which is hereby 
incorporated by reference in its entirety). The performance of a traditional linear 
hollow fiber module is compared with the Dean vortex helical hollow fiber 
module (U.S. Patent RE 37,759 to Belfort, which is hereby incorporated by 
reference in its entirety). 

Example 6 - The Aggregate Transport Model for Poly-disperse Suspensions 

[00105] As described above, the goal of the model provided by the present 
invention is to predict the performance of microfiltration of poly-disperse 
suspensions in terms of permeation flux and yield of a target species. The 
simplifying assumptions in this model are laminar flow, absence of inter-particle 
and particle-to-membrane interactions. The first step is to establish the particle 
size distribution of the suspension. Existing back-diffusion and inertial lift laws 
are then employed to calculate the hypothetical mono-disperse permeation fluxes 
for each particle size and concentration (Belfort et al., 'The Behavior of 
Suspensions and Macromolecular Solutions in Crossflow Microfiltration," J. 
Membr. ScL 96: 1-58(1 994), which is hereby incorporated by reference in its 
entirety). The lowest of these permeation fluxes is then considered the 
determining flux of the poly-disperse suspension. This permeation flux is then 
used in the back-transport laws to calculate the concentration of each species in 
the filter cake. Essentially, this is the equiUbrium concentration at the membrane 
wall that can ensure a balance between forward and back-transport of each species 
firom the membrane. The evaluated packing densities of various particles are then 
tested with respect to packing constraints that limit the cake depending on the 
particle sizes. If the packing constraints are not satisfied, the highest packing 
density is lowered and the steps executed once again. This is repeated xmtil all the 
packing constraints are satisfied. Thus, the nature of the filter cake is evaluated. 
The intOTstitial gap between the particles is estimated and steric arguments are 
used to estimate the yield of the target species (Zeman et al., ''Microfiltration and 
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Ultrafiltration Principles and Applications^' New Yoik: Marcel Deldcer, lac. 
(1996), which is hereby incorporated by reference in its entirety). If the yield of 
the target particle is between 0 and 0.94 for four diavolumes, the non-retentive 
stagnant film model is employed for the transmitted species and all the steps are 
r^eated to evaluate the corrected flux and yield. 

Example 7 - Microfiltration Characterization With Transgenic Goat Milk 

[00106] Transgenic goat milk to be used as feed su^ension was supplied 
by GTC Biotherapeutics (CHiarlton, MA). The average composition of the 
transgenic goat milk is shown in Table 4, below (Dairy Processing Handbook. 
Tetra Pak Processing Systems, AB, S-221 86, Lund Sweden, (1995), which is 
hereby incorporated by reference in its entirety). The hiunan IgG concentration in 
flie transgenic goat millc 8 g/1) was diluted with non-transgenic milk to between 
1.75 to 3 g/1. 



TABLE 4 



Composition (wt.%) 


Salient Properties 


Main Component 
Diameters 


Fat: 3.5% 


pH: 6.6-6.9 


Fat: 1 to 20 pm 


Proteins: 3.1% (80% 


Human IgG pi: 5.5 - 6.5 


Casein micelles: 


casein, rest a-lactalbu3min. 


Casein pi: 4.6 


0.3 - 0.4 pm 


lactoglobulin. 


IgG: 20mn(155 kDa) 


Immunoglobulins) 




Lactose: 4.6% 




Other whey proteins: 




15- 40 kDa 


Ash: 0.8% 






Human IgG: 2 to 3 g/1 






Total solids: 12%. 
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Example 8 - Membrane Modules 

[00107] Tubular hollow fiber membrane modules were provided by 
Millipore Corporation (Bedford, MA). Each module had six 0.1 pm mean pore- 
size poly(efher sulfone) hollow fibers. For the helical module, the fibers were 
wound in a single-wrap helix around an acrylic rod as shown in Figure 7 (U. S. 
Patent RE 37,759 to Belfort, which is hereby incorporated by reference in its 
entirety). The purpose of the rod is to stabilize the membrane helix and establish 
the desired curvature (a/r^. Table 5 summarizes the design specifications and 
operation parameters for the helical and linear hollow fiber modules. 
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TABLE 5 



Set 


Model 


■"■tn 






Linear 










L-1 


110698-1 


153 


64 


52 




1 1 n6QR-2 


153 


64 


51 




110698-4 


72 


30 


60 






72 


30 


59 


L-3 


110698-6 


44 


18.5 


86 




A X\joyo~o 


44 


18 5 


92 


HeUcal<^^'^> 










H-1 


093098-4 


141 


58.9 


38 




091198-4 


141 


58.9 


25 


H-2 


110398-2 


60 


24.9 


47 




110398-5 


60 


24.9 


49 


H-3 


111298-1 


32 


13.5 


99 




111298-2 


32 


13.5 


99 



Notes: 

a) AJs^n^) is the membrane surface area 

b) X(cm) is the active fiber length 

c) £p(hnh/kPa) is the hydraulic permeability 

d) Nominal pore diameter = 0.1 ^im 

e) Number of fibers per module, m = 6 

f) Fiber inside diameter, dj = 1.27 mm 

g) Fiber outside diameter, do =1.92 mm 

h) Rod diameter, drod for the helical modules = 6.35 mm 

i) Pitch of helical module, p = (m/ln)do = 1 .83 mm 

j) Radius of curvature of the tube in one plane, r— {drod + d^l2 = 4.14 mm 

k) Radius of curvature, = (r^ + P^lr= 4.95 mm 



Example 9 - The Microfiltration Pilot-Plant System 

[00108] Figure 8 is a flow diagram of the microlBltration pilot-plant system. 
This system consisted of two parallel flow lines, which allowed comparison 
between the linear and helical modules simultaneously for e^eriments relating to 
flux. The dual module microfiltration system consisted of a 10 L polypropylene 
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tank (Nalgene, Chicago, IL), a diaphragm pump (Mod # F 301010110, Flojet, 
Irvine, CA) and the two membrane modules. A bypass line was used to control 
the inlet flow rate to a union *T", which divided the feed flow into two parallel 
lines, one for each module. Needle valves were placed in front of the modules 
and on the downstream side to control both pressure and flow rate independently 
in each line. The inlet pressure of the modules, the transm^brane pressure 
(TMP) an the pressure drop across the modules (transcartridge pressure) as well as 
the permeate pressures were measured with digital pressure gauges (Mod # 
68920-36, PSrrronix, Tulare, CA). The retentate and permeate flow rates were 
measured separately for each fluid stream (flowmeter tube size #14 for retentate 
streams and #12 for permeate streams, Gihnont, IL). For deionized water tests 
(first phase of experimentation), a 0.2 |xm mean pore-size prefilter (Mod PSCL — 
SI, Ametek, Sheboygan, WI) was used to reduce fouling due to traces of colloids 
in the water. Two solenoid valves (Type 211, Burtek, Ingelfingen, Germany) 
located in the permeate line of each module enabled selection between normal 
operation and constant flux^ackwashing procedures. These valves were operated 
from a control xmit (GraLab Model 900, Dimco-Gray Co, Centerville, OH). Two 
peristaltic pumps (Masterflex, 7521, Cole Parmer, Chicago, IL), installed on the 
permeate lines, wcto used for either backflushing of the modules or operating at 
constant flux. For the transgenic goat milk experiments, to reduce the dead 
volume in the system, the 10 L reservoir was replaced by a 1 L graduated flask 
and the pre-filter was removed from the circuit. 

Example 10 - Characterization of Modules 

[00109] Twelve hollow fiber modules (6 linear and 6 helical; six pairs, each 

pair differing in hollow fiber membrane length) were characterized using 
deionized water with respect to flow parameters. For each module, axial pressure 
drops were varied and the corresponding axial flow rates recorded. This was first 
done wifli the permeate ports closed. Later, this was repeated with permeate ports 
opened to give a flow rate of 40 ml/min. Hydraulic permeability experiments 
were then conducted by varying the TMP from 5 kPa to 70 kPa. For all the 
hydraulic permeability tests, the axial flow was maintained at 0.5 1/min 
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corresponding to an exit Reynolds number of 1075. Two similar pairs (based on 
the length of the hollow fiber membranes) comprising of a linear and a helical 
module were selected for further experimentation with milk. 

Example 11 - Microfiltration Characterizatioii With Transgenic Goat Milk 

[00110] This phase of experiments were conducted to xmderstand and 
compare the permeation flux and protein transmission behavior of whole 
transgenic goat milk of the linear and helical modules at diffCTent Reynold's 
numbers (shear rates) and protein concentrations of milk. To study the effect of 
protein concentration, experiments were run at a fixed Reynolds number with milk 
samples correspondiag to 3, 2, 1, 1/2, and 1/3 times the normal milk 
concentration. Dilutions were achieved with addition of DI water and 
concentrations were obtained by filtering at a very low TMP (15 kPa), To 
maintain a fixed protein concentration in the bulk, both retentate and permeate 
were completely recycled to the feed reservoir. In these experiments, the TMP 
was gradually raised firom a low value to the pressure independent region of 
permeation flux. Thus, flux variation with respect to TMP could also be obtained. 
Each experiment was conducted separately for a helical and a linear module. The 
TTiillc samples were preheated to 25**C in a water batih. Samples (5 ml) were taken 
from both the retentate and the permeate and analyzed for analyzed for protein 
concentration by the Bradford (proteia) assay. 

Example 12 - Diafiltration Experiments 

[00111] Milk at concentrations corresponding to 1, 1.5, 2, 2.5, and 3 times 
the normal milk concentration were used for constant volume diafiltration 
experiments up to 5 diavolumes. For IX concentration, diafiltration was started 
after flushing the system with milk. For 2X concentration, one system volume 
was collected prior to diafiltration, while maintaining a constant reservoir level 
with milk addition. The permeate collection volumes for 1 .5X, 2.5X, and 3X 
were 1/2, 1.5, and 2 times the system volume, respectively. For tiie first protocol, 
collected permeate was recycled to the feed along with DI water in the ratio 1:1. 
These experiments were run at a fixed Reynolds number and at 90% of the 
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pressure corresponding to the pressure independent flux for the linear and helical 
modules. To study the effect of Reynolds number (shear rate), additional 
diafiltration experiments were conducted at different Reynolds numbers with the 
same milk concentration of 2X. For these experiments, the pOTneate was not 
recycled. Therefore, for these experiments concentration of the feed occurred 
during the first filtered diavolume. Samples (5 ml) were taken firom the retentate 
and permeate streams at regular intervals and analyzed for protein and IgG 
concentrations. 

Example 13 - Cleaning Protocol 

[00112] The following cleaning protocol was used. After each experiment, 
the entire system was rinsed with deionized water at an axial flow velocity of 2 
m/s for 5 minutes with the permeate ports fully opened. This was followed by 
recycUng cleaning agents Ultrasil 10 ~ detergent at 0.5 wt. % and Ultrasil 02 — 
surfactant at 0.1 wt. % at an axial velocity of 2 m/s at 45 °C for 30 minutes. The 
cleaning agents (Ultrasil 02, 10, Ecolab, St Paul, MM) were then flushed firom of 
the system for 10 minutes with deionized water. This was followed by 
sterilization with 0.1 wL % NaOCl at 40 °C for 10 minutes at 0.33 m/sec. This 
low velocity was chosen to give sufficient residence time for the bleach to act on 
the membrane modules. The membranes were stored in this dilute bleach solution 
till the next experiment. Prior to a new experiment the dilute bleach solution was 
flushed out by rinsing with deionized water for 10 minutes at 2 m/s velocity. 

Example 14 - Assay Protocols 

[00113] The Bradford assay (Bio-Rad, Hercules, CA) was used to 
determine protein concentration. Bovine lyophilized casein powder (Sigma, St. 
Louis, MO) was used as a standard and readings ware taken in disposable 5 ml 
polystyrene cuvettes (Bio-Rad, Hercules, CA). The absorbance readings with the 
spectrophotometer (Hitachi, Japan) were taken in the visible range at 595 nm 
wavelength. 

[00114] IgG assay was based on the protocol provided by GTC 
Biotherapeutics (Framingham, MA). Briefly, a protein A affinity chromatography 
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(PA ImmimoDetection™ sensor cartridge (2.1 x 30 mm) (PCTSeptive Biosystems, 
Framingham, MA) was used to obtain IgG concentrations in the various goat milV 
streams, 1 .5 ml of milk samples were pipetted into 2 ml Eppendorf centrifuge 
tubes and centrifuged at 21000g for 30 min. The milk separated into a top fat 
layer, a clear whey solution, and a casein pellet 0.75 ml of the clear whey phase 
was carefully extracted with a pipette after puncturing the fat layer. This was 
pipetted into centrifuge tubes (Spin-X tubes. Coming, NY) with 0.45 pm pore size 
cellulose acetate membranes and was centrifuged at 21000g for 15 min. The 
clarified permeate was then injected into the HPLCcolumiL For the permeate 
samples, sample preparation was unnecessary. A HPLC (Waters 510 with 
Millenniimi 2010 operating system) witiii a 486 UV detector and U6K sample 
injector ware used (Waters Corp., Milford, MA). The loading buffer was 1 0 mM 
phosphate bxiffer and 150 mM NaCl at pH 7.20 ± 0.05, while the elution buffer 
was 12 mM HCl with 150 wM NaCl. The pump flow rate was set at 2 ml/min., 
and the detector wave length at 280 nm. The injection volume was 10 \iL for milk 
and 20 to 40 for permeate samples. A calibration graph was constmcted by 
injecting different dilutions of IgG fusion protein (GTC Biotherapeutics, 
Framingham, MA). Loading buffer was passed through tibie colunm for 1 0 
minutes followed by sample injection and loading buffer again for 5 minutes. 
After tins, elution buffer was run for 10 minutes. A clean peak corresponding to 
IgG fusion protein was detected at around 6.5 minutes into the elution phase. 
Area obtained by peak integration was compared with the calibration graph to 
obtain flie IgG concentration of the sample after dividing by the sample volume. 
Care was taken to ensure that all readings were within the range of the calibration 
graph. This was done by adjusting the sample injection quantities. 
[00115] Fat content was measured by the Gerber method which is approved 

for use by dairies in USA. Eleven ml of preheated milk sample (37^C) was added 
to 10 ml of sulfuric acid in a butyrometCT. 1 ml of amyl alcohol was added, and 
the butyrometer was capped with a special stopper. Shaking the butyrometer 
CTisured digestion of the proteins by sulfuric add. The butyrometer was then 
inverted and centrifiiged for 6 ndnutes at 350g. After this, the butyrometer was 
immersed in water bath at 65°C for 5 min. The fat appeared as a clear liquid, and 
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ihc quantity was read out as a volume percentage in tibie graduated section of the 
butyrometer. 

Example 15 - Model Predictions 

[00116] Model predicted pressure-independent pemieation fluxes during 
microfiltration of whole transgenic goat milk at 298 K, with a 6-fiber hollow fiber 
module of length 300 mm, internal diameter 1.27 mm^ and pore diameter of 100 
nm, were plotted against mean axial shear rate for mono-disperse suspensions of 
IgG, casein micelles, and fat globtdes, shown in Figure 9A. The permeation 
fluxes were based on back-transport of these particles from the cake at the wall. 
The mi'niTnum of these fluxes was used as the first estimate for the permeation 
flux of whole transgenic goat milk which is a poly-disperse suspension 
comprising these three particles as the dominant specie. This value was further 
refined by successive iterations in the model to the final value of the pressure- 
independent poly-disperse permeation flux for whole transgenic goat milk. This is 
plotted against mean axial shear rate, shown in Figure 9B. From Figure 9A, it is 
seen that for shear rates below 40,000 s"\ the IgG molecules exhibit higher back- 
transport than the casein micelles, a casein cake layer is expected to limit the 
pressure-independent permeation flux. A high value of IgG transmission is 
expected in the casein cake regime as the interstitial gap between the casein 
micelle spheres will be > 150 nm which is much higher than 20 tun the diameter 
of the IgG molecules. For shear rates higher than 40,000 s"\ the model predicts 
that the cake layer will be mainly composed of IgG. This will form a tighter layer 
resulting in reduced IgG transport. 

Example 16 - Microfiltration Characterization With DI Water 

[00117] Experiments were perfomied with DI water for 12 modules as p&r 
Table 5 to determine the hydraulic permeability and fiiction factor. It is seen in 
Figure 10, that the friction factor followed the well known Hagen Poiseuille 
equation for the linear modules and the Mishra-Gupta equation for helical 
modules. The intermediate length modules, L-2 for linear and H-2 for helical 
were selected for further experiments with transgenic whole goat milk. 
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Example 17 - MicroflUitration Characterization With Whole Goat Milk 

[00118] The data in Figure 1 1 and Figure 12 indicate that the classical gel- 
polarization model appears to be valid for milk filtration. This was further 
confirmed by a plot of flux versus the logarithm of the protein feed concentration 
which yielded straight lines for both modules, shown in Figure 13, The data 
indicate a significantly higher mass transfer coefficient for the helical module as 
compared to tiie linear module (26.5 versus 18.2 Imh) as expected (Luque et al., 
"A New Coiled Hollow-Fiber Module Design for Enhanced Microfiltration 
Performance in Biotechnology," Biotechnol Bioeng, 65:247-257 (1999); Belfort 
et al., *The Behavior of Suspensions and Macromolecular Solutions in Crossflow 
Microfiltration," J. Membr. ScL 96:1-58 (1994); U. S. Patent RE 37,759 to 
Belfort, which are hereby incorporated by reference ia their entirety). The plot is 
made according to the gel polarization model J = k Ln(Cw/Cb). Based on the x- 
intercept, Cw (helical) == 448 g/1 and Cw (linear) = 443 g/1. These values compare 
well with values reported in literature of 280 g/1 for whey protein filtration 
(Zeman et al., ^^Microfiltration and Ultrafiltration Principles and Applications^^ 
New York: Marcel Dekker, Inc. (1996), which is hereby incorporated by reference 
in its entirety). The transmembrane pressures corresponding to the pressure- 
independent flux, were about 60 kPa for the linear and 70 kPa for the helical 
modules. 

Example 18 - Diafiltration Experiments 

[00119] A series of diafiltration experiments were conducted to study the 
IgG transmission and permeation flux behavior with the passage of diafiltration. 
Experiments at a fixed Reynolds number of 1400 and diflEerent starting milk 
concentrations indicated similar flux trends for both modules, shown in Figure 14, 
For each case, the heUcal module gave a 55 to 70 % higher flux compared with 
the linear modules. For lower concentrations, IX and L5X, tiie flux increased 
with diavolumes. This behavior is e^qplained by two competing phenomena. 
There was a progressive reduction of viscosity in the bulk solution due to the 
reduction in total soUds content especially the sugars, salts, and low molecular 
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weight moieties. This enhanced the back-transport of casein as well as IgG, the 
putative foulants according to an earlier hypothesis, and a higher forward 
convective transport of solvent could be supported without cake build-up. For 2X 
and 2.5X concentrations the flux remained steady through most of the run witii a 
drop at the end. This was probably due to the cancellation of the deleterious effect 
of cake build-up and the beneficial effect of reduction in suspension viscosity. 
For higher diavolumes, cake build-up overwhelmed the effect of viscosity 
reduction and a decreasing trmd was obtained at the end of the run. This is 
expected, because viscosity reduction tapered off with filtration as the total solids 
firaction became constant due to low sieving at the end of diafiltration. For 3X 
concentration, there was a continual reduction of permeation flux with 
diafiltration during the run. This is explained by the overriding influence of cake 
build-up due to poor back-transport as well as high viscosity of the suspension. It 
was also foxmd that the IgG and whey protein transmission was lower for the 
helical modxile, shown Figure 15 and Figure 16. 

[00120] hi the second phase of diafiltration, shown in Figure 1 7, it was 

found that for each case there is an increase in permeation flux with Reynolds 
number. As indicated earli^, in these experiments, conducted with a milk 
concentration of 2X, the permeate was not recycled. Hence, the first diavolmne 
comprised the concentrating phase and as expected, a significant drop in 
permeation flux was observed during this phase. The other interesting point is 
that for hi^ Reynolds number flow the permeation flux rose afl:er the 
concentration phase for the helical module, hi all other cases, the flux dropped. 
This is explained by the fact that at the same Reynolds number the helical module 
has a wall shear rate of almost twice that of the linear module. At these high shear 
values, the cake build up was slow or non-existent and the increase in permeation 
flux was due to the reduction in suspension viscosity. Figure 1 8 indicates that the 
helical module had a permeation flux of 55 to 70 % higher than the linear module 
for the same Reynolds number. 
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Example 19 - Fit of Experimental Data and Predictions 

[00121] A very good fit (r^ = 0.92) between the model and experiments for 
a wide range of operating conditions (near the pressm'e-independent flux regime) 
with variations in Tnilk concentration, temperature, and Reynolds numbers was 
obtained, as seen in Figure 19. For the helical module, in all cases except at a 
Reynolds number of 640, the shear rate is high^ dian 40,000 sec*\ indicating that 
the IgG molecules will build up the cake instead of casein micelles. This explains 
the lower protein transmission but high flux while opCTating under high wall shear 
with the helical module. Operating at high shear rates, but at a permeation flux 
which is significantly lower than the value predicted by the model is 
recommended. In addition, to prevent increase in permeation flux during 
diafiltration at constant pressure, filtration at constant permeation flux rather than 
constant transmembrane pressure is suggested. Under these conditions a cake wiU 
likely not form completely on the membrane and good protein transmission of the 
target molecule, IgG, during nnlk naicrofiltration in the diafiltration mode will be 
observed. 

[00122] This concept has been addressed by others (Field et al., "Critical 
Flux Concept for Microfiltration Fouling," J. Membr. Set. 100:259-272 (1994); 
Ould- Dris et al., ''EflFect of Cake Thickness and Particle Polydispersity on 
Prediction of Permeate Flux in Microfiltration of Particulate Suspensions by a 
Hydrodynamic Diffusion Model," J, Membr. Set 164:21 1-227 (2000); Gesan- 
Guiziou et al., "Critical StabiUty Conditions in Skimmed Milk Crossflow 
Microfiltration: Impact on Operating Modes," Zraft 80:129-140 (2000), which are 
hereby incorporated by reference in their entirety). The model allows the 
establishment, a priori, of a safe operating permeation flux. This can be 
generalized to microfiltration of other well characterized poly-disperse 
suspensions to arrive at the safe permeation flux. This also underscores the 
concept of constant permeation flux operation as preferable to constant pressure 
microfiltration. This is because with a constant pressure approach, it is likely that 
the safe value of permeation flux, above which a cake layer forms completely, 
will be exceeded during the initial stages of tiie run and during the course of 
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diafiltration due to the reduction of protein volume firaction and bulk viscosity due 
to lower volume firaction of solids in the bulk feed. 

[00123] A series of experimental results with whole transgenic goat m ilk 
under various operating conditions of temperature, wall shear rate, and milk 
concentrations were presented herein. These results validate the predictive 
aggregate transport model for microflltration of combined macromolecular 
solutions and poly-disperse suspensions of the present invention. Milk, being an 
extremely complicated poly-disperse suspension, provided a suitable test 
suspension for the method disclosed herein. Further to the recommendations 
herein, salient features and benefits that can be derived from this model. 
Crossflow microfiltration performance of different poly-disperse suspensions can 
be predicted a priori. Currentiy, in the absence of a general theory, numerous 
experiments need to be p^formed in the laboratory and bench scales for each feed 
type. Therefore, using the model, there is a tremendous potential for saving time 
and labor. The model can be used with different geometries. The model can 
determine the nature of the filter cake. The model, being theoretical, can be used 
for scale-up and scale-down for industrial or laboratory applications. This is a 
major pitfall of empirical methods which are vaUd only for the operating region 
and scale. The computerized version of the model can be interfaced with other 
software packages for optimizing diafiltration for the optimum plant operation. 
[00124] Although preferred embodiments have been depicted and described 
in detail herein, it will be apparent to those skilled in the relevant art that various 
modifications, additions, substitutions, and the like can be made without departing 
from the spirit of the invention and these are therefore considered to be within the 
scope of the invention as defined in the claims which follow. 



